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1.  General  Remarks* 

During  the  period  November  1,  1981  to  October  31,  1982,  the 
Principal  Investigator,  in  cooperation  with  several  research 
assistants,  carried  out  a  orogram  of  mathematical  research  in 


the  general  area  of  control  theory  of  partial  differential  equa¬ 


tions.  The  program  involved  two  distinct  phases;  an  effort  aimed 


specifically  at  the  development  and  improvement  of  control  strat¬ 
egies  in  connection  with  the  wing  flutter  problem  and  a  more 
general  program  in  the  area  of  distributed  parameter  control  prob¬ 
lems  of  hyperbolic  type. 

This  work  resulted  in  two  scientific  papers  which  form  the 
greater  part  of  this  report-  The  first  of  these,  "Some  Remarks 
on  the  Current  Status  of  the  Control  Theory  of  Single  Space  Dim¬ 
ension  Hyperbolic  Systems"  was  presented  at  the  NASA  JPL  Symposium 
on  Control  and  Stabilization  of  Large  Space  Structures,  Pasadena, 
CJK,  July,  1982.  The  second,  "Admissible  Input  Elements  for  Systems 
in  Hilbert  Space  and  a  Carleson  Measure  Criterion",  by  £.  P.  Ho 
and  the  Principal  Investigator,  is  a  paper  which  largely  resulted 
from  Dr.  Ho's  thesis  work,  also  supported  by  this  grant,  in  part. 


In  addition  to  Dr.  Ho,  who  is  now  with  the  University  of  Iowa, 
the  Principal  Investigator  was  assisted  by  R-  G.  Teglas,  H.  M. 
Baron,  and  R..  Rebarber.. 
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2. 

2*.  Travel  Supported  by  the  Grant. 

Grant  funds  were  used  to  support  travel  by  the  Principal 
Investigator  and  one  Research  Assistant,  H.  U.  Baron. 

The  Principal  Investigator  travelled  to  Pasadena,  California, 
to  attend  and  present  a  paper  at  the  NASA  Jet  Propulsion  Laboratory 
Symposium  on  Control  and  Stabilization  of  Large  Space  Structures. 
The  paper  presented  concerned  the  current  status  of  the  control 
theory  of  hyperbolic  partial  differential  equations  with  particular 
emphasis  on  observers  and  canonical  structure. 

The  Principal  Investigator  also  took  part  in  the  30th  Anni¬ 
versary  meeting  of  the  Society  for  Industrial  and  Applied  Mathe¬ 
matics  (SIAM)  in  Palo  Alto,  California,  July,  1982.  A  paper  out¬ 
lining  the  treatment  of  control  problems  associated  with  infinite 
dimensional  linear  systems  by  means  of  methods  from  the  theory  of 
analytic  functions  was  presented  at  this  meeting  by  invitation  of 
the  organizing  committee*.  Ms..  Baron  also  attended  this  meeting 
and  presented  a  paper  on  control  canonical  forms  for  systems  gov¬ 
erned  by  various  types  of  partial  differential  equations. 


3-  Technical  Appendix 


This  appendix  consists  of  two-  papers  whose  preparation  was 

supported  in  part  by  the  grant.  These  papers  ares 

"Some  Remarks  on  the  Current  Status  of  the  Control  Theory 
of  Single  Space  Dimension  Hyperbolic  Systems" 

and 

"Admissible  Input  Elements  for  Systems  in  Hilbert  Space 
and  a  Carleson  Measure  Criterion", 

the  latter  paper  being  jointly  authored  by  L.  P.  Ho  and  the 

Principal  Investigator. 


SOME  REMARKS  ON  THE  CURRENT  STATUS  OF  THE  CONTROL  THEORY 
OF  SINGLE  SPACE  DIMENSION  HYPERBOLIC  SYSTEMS 

D.  L.  Russell’* ** 

University  of  Wisconsin  Mathematics  Department 

ABSTRACT 

We  review  various  aspects  of  the  control  theory  of  hyperbolic  systems, 
including  controllability,  stabilization,  control  canonical  form  theory,  etc.  To 
allow  a  unified  and  not  excessively  technical  treatment,  we  restrict  attention 
to  the  case  of  a  single  space  variable;  the  multi -dimensional  case  is  treated  in 
our  more  extensive  review  [  36]  .  The  paper  concludes  with  a  short  discussion 
of  the  newly  developed  procedure  of  canonical  augmentation. 

SOME  ASPECTS  OF  THE  CONTROL  THEORY  OF  THE  WAVE  EQUATION 

AND  RELATED  SYSTEMS 

The  systematic  study  of  control  systems  governed  by  partial  differential 
equations,  a  special,  but  exceptionally  important,  subcategory  of  distributed 
parameter  systems  began  in  the  early  1960's  with  the  work  of  the  Soviet 
scientists  A.  G.  Butkovskii  [  3  ]  ,  [4],  Yu.  V.  Egorov  [11]  and  others.  These 
works  were  primarily  concerned  with  the  extension  of  Pontyagin*  s  Maximum 
Principle  [26]  to  certain  classes  of  processes  which  could  not  be  satisfactorily 
modelled  by  finite  dimensional  mathematical  systems.  Controllability  questions 
were  raised  but  were  usually  subsidiary  to  questions  of  optimality.  One  of  the 
first  systematic  controllability  studies,  in  connection  with  the  heat  equation, 
was  presented  by  Gal'chuk  in  [14]  .  One  of  the  most  important  of  the  early 
American  contributions  to  the  subject  was  the  1963  thesis  of  Fattorini  [13]  , 
which  also  treated  parabolic  systems  and  was  one  of  the  first  works  to 
recognize  the  strong  relationship  between  distributed  parameter  control  studies 


*  Supported  in  part  by  the  Air  Force  Office  of  Scientific  Research  under  Grant 
AFOSR  79-0018. 

**Dept.  of  Mathematics,  University  of  Wisconsin,  Madison,  WI  53706.  Also 
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and  classical  results  in  analytic  function  theory. 

The  author*  s  own  interest  in  distributed  parameter  control  theory  arose 
out  of  consulting  experience  with  Honeywell,  Inc. ,  and  NASA,  starting  around 
1965  or  1966.  In  developing  the  Saturn  launch  vehicle  for  the  Apollo  program, 
NASA  has  encountered  the  problem  of  transverse  vibrations  of  the  booster 
structure  and  interaction  of  those  vibrations  with  liquid  sloshing  modes  in  the 
immense  Saturn  fuel  tanks.  While  the  eventual  treatment  of  that  problem  was 
based  on  finite  modal  approximations,  the  problem  stimulated  a  great  deal  of 
research  aimed  at  an  understanding  of  the  control  of  vibrations  in  various 
distributed  parameter  settings. 

First  looking  at  this  problem,  under  Honeywell -NASA  auspices,  we 
thought  of  modelling  the  booster  structure  as  an  "Euler**  beam,  the  displacement 
w(x,  t),  which  we  may  take  to  be  scalar  here,  satisfying 

p(x)  ~  +  -^-  ( El (x)  )  a  0  (1.1) 

at*  ax*  ax* 

along  with  appropriate  boundary  conditions  including  the  control  inputs,  at  the 
longitudinal  extremities  x  =  0,  x  =  L  .  We  got  nowhere  with  our  study  of 
this  problem  initially  because  the  equation  (1. 1)  is  not  particularly  well 
understood  from  the  mathematical  standpoint.  There  seemed  to  be  no  "handles" 
to  grasp.  It  would  not  be  until  the  1969  thesis  of  Quinn  [  27]  that  we  would 
understand  how  this  system  works  and  that  it  is,  in  fact,  controllable  in  a 
rather  strong  sense. 

We  knew  about  the  control  theory  of  ordinary  differential  equations  from 


various  papers  and  from  notes  and  lectures  which  would  later  be  incorporated 
into  the  1967  treatise  on  control  theory  by  Lee  and  Markus  [20]  .  We  also 
knew  that  hyperbolic  partial  differential  equations  in  two  independent  variables 
reduce  to  ordinary  differential  equations  satisfied  along  the  characteristics.  It 
was  natural,  therefore,  to  look  for  hyperbolic  models  which  might  fit  our  purpose. 
Such  was  provided  by  the  Timoshenko  beam  equations 

<l-Z) 


(1.3) 


which  may  be  viewed  as  two  coupled  wave  equations.  By  "wave  equations" 
here,  we  mean  the  equation 


-  -irUM-lS-)  =  o. 


(1.4) 


All  coefficient  functions  shown  in  (1. 2),  (1.  3),  (1.  4)  are  positive  on  OsxsL, 

It  may  be  verified  that  (1. 2),  (1. 3)  and  (1. 4)  are  hyperbolic  in  the  sense 
described  in  £  8]  ,  [25]  ,  e.g. .  Since  (1.4)  is  conceptually  simpler,  it  was 
studied  first,  with  accompanying  boundary  conditions 

z(0,t)  -  0  (1.5) 


f(I*t)  =  u(t). 


(1.6) 


the  latter  incorporating  the  control  force  u  (t)  . 

While  the  practical  goal  in  mind  was  appropriate  form  of  stabilization, 
we  knew  that  in  the  case  of  finite  dimensional  systems 

x  =  Ax  +  Bu 

an  affirmative  resolution  of  the  controllability  problem,  steering  from  a  given 
x(0)  a  Xg  to  a  given  x(T )  *  x^  implied  the  property  of  stabilization;  hence 
we  felt  Justified  in  first  looking  at  the  state  to  state  controllability  problem  for 
(1. 4),  (1. 5),  (1.  6).  The  "energy"  form  for  (1. 4)  is 


(1. 4),  (1. 5),  (l.  6).  The  "energy"  form  for  (1. 4)  is 

e(*-  ir  > ■  ?  X  r«  <  tI  <*•  *  »2  +  s(x>  <  -fe-  <*■ l>  ,2<J* 

Given  initial  and  terminal  states 

z(x,0)  =  zn(x),  4l(x,  0)  =  yn(x) 


(1.7) 


(1.8) 


Z  (  X,  T  )  =  Z|(X  ),  (X,T)  a  yx(x)  (1.9) 

of  finite  energy,  i.e.  e(zQ,  v0 )  <  •  ,  C(Zp  vj)  <  •  ,  we  asked  if  there 
exists  u  e  L^[0,  T]  for  which  the  solution  of  (1.4),  (1.5),  (1.6)  correspond¬ 
ing  to  the  initial  state  (1.8)  assumes  the  desired  terminal  state  (1.9)  at  time 
t  =  T  .  The  answer,  a  qualified  "yes",  came  from  two  different  approaches  to 
the  problem.  The  relationship  between  these  two  approaches  has,  over  the 
years,  grown  ever  more  fundamental  and  has  led  to  a  great  many  very  interesting 
developments.  See  [34]  and  [45]  in  particular. 

The  first  method  explored  was,  as  we  have  already  indicated,  the  method 
of  characteristics.  If  we  let 

n- io) 

and  consider  families  X+,  X"  of  "characteristic"  curves  satisfying 


(1. 10) 


-fp-  +  C  (x)  =  0,  -  c  (x)  =  0 ,  (1. 11) 

respectively,  and  then  set 

v+(x,  t )  =  -|^*  (x,  t )  +  c  (x)  (x,  t ) , 

v“(x,t)  =  |-  (x.t)  -  c(x)  ~  (x.t)  , 

we  see  readily  that  on  X+=  {(x+(t),  t)>  ,  X“  =  {(x-  (t),  t)>,  respectively, 
we  have 

—  v+(x  +  (t),  t)  =  c'(x  +  (t))  Y  (x  + (*)>*?  ~  v  (x  +  (t),t) 

2 

v“(x-(t),  t)  s  c'(x-(t))-  tx -(*)>*)  ~  y"(.x~  t?  (1.13) 

2 


t  “  t  >  2Tl 


Fig.  1. 1:  The  Method  of  Characteristics 


r^rt avtTv lv .»; -  v.vv  «/  wr  v  -.■*.  \^. 


Because  these  differential  equations  are  satisfied  on  different  families  of 
characteristics,  the  coupling  between  them  is  more  complicated  than  for  the 
usual  system  of  ordinary  differential  equations.  Nevertheless  there  is  a  method 
of  successive  approximations,  described  in  [30],  [36],  which  enables  solution 
of  these  equations  in  certain  regions  provided  with  appropriate  boundary  data. 
Such  a  region  is  the  roughly  triangular  domain  Aq  shown  in  Fig.  1. 1.  bounded 
by  t  *  0,  x  =  0  and  the  characteristic  X+(L,  0),  of  the  first  family 
described  by  (1.11),  passing  through  the  point  (L,  0).  Together  with  the 
boundary  data  provided  by  (1. 5)  and  (1.  8),  it  may  be  seen  that  the  differential 
equations  (1.12),  (1.13)  determine  v+  and  v“,  and  hence  z(x,t),  through¬ 
out  the  domain  Aq  .  Similarly,  these  equations  together  with  the  data 
provided  by  (1.5)  and  (1.9)  determine  z(x,  t)  in  the  domain  Aj  bounded 
by  x  *  0,  t  =  T  and  the  characteristic  curve  X~(L,T),  described  by  the 
second  equation  in  (1. 11)  and  passing  through  the  point  ( L,  T )  .  Thus  the 
initial  and  terminal  states,  described  by  (1. 8),  (1.9).  together  with  the 
boundary  condition  (1.5)  determine  z(x,  t)  in  both  AQ  and  a^. 


Whether  aq  and  A^  are  disjoint,  or  have  a  region. 


Qr 


of  over¬ 


lap  ,  depends  on  the  time  T  allotted  for  control.  The  time  required  for  the 

curve  X  (0,  L)  to  pass  from  X  *  L  to  x  *  0  is 

_L 


Ti  = 


l  (U4) 

and  this  is  also  the  time  required  for  X~(L,T)  to  pass  from  x  =  0  to  x  =  L. 
We  summarize  the  control  situation,  depending  on  the  relationship  between  T 
and  Tj . 

Case  T  <  2T  .  Here  AQ  and  A^  overlap  and  the  determinations  of 
z(x,t)  in  the  overlap  region  nQ=  AQ  U  A^  provided  by  (1.8)  and  (1.9)  need 
not  and,  in  general,  will  not  agree.  There  can,  in  such  cases  of  disagreement, 
be  no  solution  of  (1.4),  equivalently  (1.12),  (1.13),  in  the  region 
Rj,  s  {(x,  t)|  0  ^  x  as  L,  OstsT}.  The  control  function  u(t),  shown  in 
(1.  6),'  never  enters  the  picture  because  it  cannot  affect  the  solution  of  (1.  4)  in 
Aq  or  Aj  if  (1.8),  (1.9)  are  satisfied  at  t  =  0,  t  =  T,  respectively. 

Case  T  =  2T  .  Here  the  two  "domains  of  determinacy",  aq  and  a^, 
just  fail  to  overlap;  their  boundaries  have  exactly  one  point  in  common,  t  =  T^, 
x  =  0  .  The  initial  and  terminal  conditions  (1.8)  and  (1.9)  determine  z(x,  t) 
in  Aq  and  A^ ,  respectively.  Another  process  of  integration  of  the  coupled 


differential  equations  (1.12)  and  (1.13)  permits  unique  extension  of  z(x,  t), 
equivalently  v+(x,t),  v“(x,  t),  into  the  domain  .  The  control  steering 
(1.  8)  to  (1. 9)  is  then  uniquely  determined  from  this  extension  and  (1.  6). 

The  determinations  of  z(x,  t)  in  and  ^  may  fail  to  match 

smoothly  at  the  point  p  :  x  »  0,  t  =  T,  .  This  results  in  discontinuities  of 

+  +  .  *  . 
v  along  X  ( L,  0 )  and  of  v  along  X  ( L,  T  )  in  general. . 

Case  T  >  2T  .  The  only  difference  between  this  case  and  the  case 

T  =  ZTj  lies  in  the  line  segment  £  :  x=0,  T^  <  t  <  T  -  T^  which  replaces 

the  point  p  of  the  case  T  =  2Tj  .  Extension  of  z(x,t)  from  U  A1 
into  Cl  cannot  be  carried  out  until  the  boundary  condition  (1. 5),  which 
yields  9z/3t  ( 0,  t )  =  0  ,  is  augmented  by  arbitrary  data 

|f(0,t)  =  C(t)  ,  (0,t)  e  £  .  (1.15) 

Once  this  is  done,  extension  of  z(x,t)  into  Cl  proceeds  much  as  before. 
(See  [30]  ,  [36]  for  details  of  the  extension  process.)  The  arbitrary 
function  r  (t )  can  be  designed  so  as  to  eliminate  discontinuities  of  the 

^  -L 

solution  along  X  (L,  0)  and  X  (L,T),  to  satisfy  some  criterion  of 
optimality  ( see  [  30]  e.  g.  )  or  to  fulfill  any  other  appropriate  design  objective. 

If  the  partial  differential  equation  (1.  4)  is  combined  with  boundary 
conditions  different  from  (1.  5),  (1.  6),  but  still  admissible  for  (1.  4),  the 

cases  T  <  2Tj,  T  >  2T^  remain  as  above.  The  rather  delicate  situation 
at  T  =  2Tj  depends  on  the  specific  form  of  the  boundary  conditions.  For 
example,  the  boundary  conditions 

z(0,t)  *  0  ,  z  ( L,  t )  *  u(t) 

lead,  in  case  T  =  2T^  to  a  situation  where  the  desired  control  is  not  unique; 
it  has  the  form 

u  ( t )  =  u(t)  +  yu(t) 

where  u(t)  is  a  non-zero  control  steering  the  zero  initial  state  into  the  zero 
final  state  and  y  is  an  arbitrary  constant.  By  contrast,  the  boundary 
conditions 

ff(0,t)  =  0.  |(U)Su{t)  (1.16) 

lead,  in  case  T  =  2T^,  to  a  situation  where  the  desired  control  u(t)  does 
not,  in  general,  exist.  (See  [31],  [37]  for  more  details. ) 

The  analysis  of  more  complicated  systems  of  hyperbolic  equations,  such 


\ 


as  the  Timoshenko  system  (1.2),  (1.3),  is  in  general  rather  compl] 
there  are  some  special  cases,  including  appropriate  boundary  condi 
which  the  analysis  is  fairly  simple.  In  [  M]  a  discussion  is  giv 
analysis  of  the  free  boundary  case 

||  (0,t)  =  0,  iMO.t)  -  U  (0,t)  =  0, 


3x 


J|  (L,  t)  =  Ul(  t)  ,  t|i(L,t)  -|J  (  L,  t )  =  u2(t)  . 

It  may  be  shown  that  all  cross-coupling  is  of  low  order  and  the  pro! 
essentially  equivalent  to  two  problems  (1. 4)  with  boundary  condit: 
Two  critical  times  are  involved.  With  (cf.  (1.14)) 

c  fxl  -  \z  t  =  f  - 

cl(x)  "  (  I  p(x)  '  '  A1  cr(x)  ’ 


c2<*>  =  T2=J[LdX 


c,(x) 


it  may  be  shown  that  finite  energy  states  are  controllable  if  and  oh 

T  >  2  max  {Tp  T2>  . 

The  essential  details  of  the  analysis  are  given  in  [30]  and  are  qi 
to  what  we  have  briefly  outlined  here  for  (1. 4),  (1. 5),  (1.  6)  . 

It  is  immediately  clear  that  the  method  of  characteristics  is 
adapted  to  controls  u(t)  acting  at  a  point,  as  in  (1.6).  This 
because  the  control  determination  occurs  at  the  very  last  stage  of  t 
after  the  controlled  solution  has  been  computed.  If  the  control 
scalar,  acts  on  the  system  through  a  "control  distribution  function 
in  (cf.  (1. 4)) 

r<*>  Tt--fe<s<x>l§>  =  B(x)u(t>- 

9t 

homogeneous  boundary  conditions  (cf.  (1.5),  (1.6)) 

•(°-*>  *  H11'1’  -  0 

applying  at  the  boundaries,  we  face  what  appears  at  first  glance  tc 
different  situation  than  what  obtains  in  (1.  4),  (1.  5),  (1.  6),  for  ev« 
equations  corresponding  to  (1.12),  (1.13)  will  involve  the  unknown 
in  this  situation;  one  cannot  proceed  by  filling  out  z  ( x,  t )  in  si 
domains  as  before;  a  completely  different  approach  is  required.  ! 
approach  can  be  found  in  the  study  of  moment  problems  -  a  techniq 


by  several  authors  (see  [3]  ,  [12]  ,  [15]  ,  [14]  ,  [23]  ).  The  technique  has 
the  advantage,  from  the  point  of  view  of  approximation  of  being  intimately 
connected  with  the  modal  representation  of  the  system  based  on  the  natural 
modes  of  vibration,  or  eigenfunctions  of  the  operator  -r(x)“*(3/3x)(s(x)(3z/3x)). 

It  is  known  (  see  [1]  ,  [7]  )  that  the  operator 


Lz  =  -r(x)"1  (s(x)  ) 


3x 


with  boundary  conditions  conformable  with  (1.  22)  has  eigenvalues 

r2 


ir 


+  e, 


k  =  1,  2,  3, 


(1.  23) 


(1.  24) 


where  the  are  uniformly  bounded  and  is  related  to  c(x)  by  (1. 14). 
The  corresponding  eigenfunctions,  ?k(x) »  k  =  1,  2, 3,  •  •  • ,  form  an  ortho- 
normal  basis  for  L2  [0,  L]  (which  consists  of  the  same  functions  as 


L2[0,L] 


L‘[O.L] 

but  has  the  inner  product 
rL 

(9,  4»)r  55  /  r (x) ? (x) <]*  (x)  dx )  . 
r  0 


(1.  25) 


Every  finite  energy  solution  z(x,t)  of  (1.21),  (1.22),  i.  e.  every  solution 

for  which  the  integral  (1.  7)  Is  bounded  for  all  t,  can  be  expanded  In  the  form 

0» 

Z(X,t)  =  Y,  2k(t)  **(*) 

k=  1 

where,  if  we  assume  the  control  distribution  function  f  (x)  has  the  expression 

aa 

g(x)  =  Y  Sk^k*** 


k=l 


2 

convergent  in  Lr[0,  L]  ,  the  z^(t)  satisfy 
zk  +  *Tc  zk  =  gk  u(t)»  k*  1,2, 


(1.  26) 


Letting 


“k  =  Nc  ' 


k  =  1,  2,  3, 


and  using  the  transformation 

l/iuk  -  l/iwk 


/2k\  fl'l‘ 

\\)  \  1 


one  arrives  at  the  system 


iw, 


(1.27) 


u(t)  .  (1.  28) 


It  may  be  seen  that  finite  energy  states  are  those  for  which 

Z  I<*ic>2  + W2]  <  - 

k  =  1 

and  this  becomes,  in  terms  of  Tk,  £k  , 

Z  tiiki2  +  icki2)  <  -  • 

k=l 

Integrating  (1. 28),  we  have,  for  T  >  0  , 

iwvT  gv  -T  iu^CT-t) 

nk(T)-e  *  nk(0)  =  -f-  /  e  u(t)dt 


(1.  29) 


u  (t)dt 


Ckft)- 


iwt,T 


gv  ^  -  iukfi-t) 

ck<°)  =  -r-^e  u(t,dt 


Assuming  the  controllability  condition 

gk  /  0,  k  =  l,  2,  3,  •  •  • 

we  see  that  the  problem  of  steering  between  the  given  states  at  times  0  and  T 
reduces  to  the  moment  problem 


f(s)ds  = 


-T  -iw-s  \ 

J  a  *  f(s)ds  =  -*-, 
0  ''k 


k  «  1.  2,  3, 


k  «  1,  2,  3. 


(1.30) 


(1.31) 


(1.32) 


where  s  =  T-t ,  f(s)  *  u  (T  -  s),  and 

1  T  T 

ak  *  2(nkfr)  "  e  ^  V0,)'  Pk  =  2Uk^>  -  e  k  Ck(0))  (L32) 

are  square  summable. 

To  solve  the  moment  problem  we  resort  to  the  theory  of  nonharmonic 
Fourier  series  as  developed  by  Paley  and  Wiener  [  24]  ,  Levinson  [  21]  , 
Schwartz  [42]  and  many  others.  (An  excellent  expository  treatment  [49]  by 
R.  Young  has  recently  appeared. )  The  following  is  known;  the  three  cases 
being  divided  in  a  manner  conformable  with  the  three  cases  discussed  earlier. 

lui,t  -i«kt 

Case  T  <  ZTj .  The  functions  e  ,  e  ,  k  *  1, 2, 3,  •  •  • ,  are 
linearly  dependent  in  L2[0,T]  in  a  rather  strong  sense.  Any  one  of  these 
functions,  indeed,  any  finite  number  of  them,  lie  in  the  closed  span  of  the 
remaining  functions  (which,  in  fact,  is  equal  to  the  whole  space  L  [0,T]  ). 
As  a  result  the  moment  problem  (1.30),  (1.31)  cannot,  in  general,  be  solved. 


I 


a 


>1 


si 


Suppose,  e.g. ,  all  but  finitely  many  of  the  3^,  say  k  =  K+l,  K  +  Z,  •••, 

wens  equal  to  zero,  while  some  of  the  ,  3^,  k  =  1, 2,  •  •  * ,  K,  are  non-zero 
The  linear  dependence  just  referred  to  shows  that  such  a  problem  can  have  no 


solution;  the  equations 


jr  w  s  _ 

J  e  *  f (s)ds  =  /  e  *  f(s)ds  =  0,  k  >  K 
0  0 

imply  that  the  same  equations  must  hold  for  k  s  K  . 


.T  iui,s 


it*>  US 


■iu^s 

e  ,  k  =  1,2,3, 


Case  T  =  Tj,  Here  the  functions  e  ,  e  ,  k  =  1, 2, 3,  •  • 

form  a  Rlesz  basis  for  L2[0,  2Tj]  .  Every  function  he  [0,  ZTj]  has  the 


unique  convergent  expansion 


iu^S 


-Iuj.s 


h(s)  =  £  thke  +  h-ke  1 
k=  1 

and  there  are  positive  numbers  c,  C,  such  that 


(1.33) 


c*zyhU2,  a  Y  <|h,l2  +  |h  k|2aC2|h[z 

i^o.zr]  "j  *  -k'  l2[0,zt] 


(1.34) 


Further,  there  is  a  unique  dual  basis  of  biorthogonal  elements  p^,  P  ^  e 


L  [  0,  2T,  ]  such  that 

1  -yn 


.“1  Ws  ks  +1,  fl,  ... 

/  «  P.Wds  *  Jl  *  , 

0  £  K,z  l  s  ±  1,  ±  2,  . . . 


*kl  = 


ri.  k  =  i 

l  0,  k  1 


(1.  35) 


which  engenders  expansions  similar  to  (1.33),  the  roles  of  c,  C  in  the 
inequalities  parallelling  (1.  34)  being  reversed.  The  formal  solution  of  (1.  30), 
(1.  31)  is  then  uniquely  given  by 


«■>.-  I  •5rC“kpk(s,  +  \p-kts)]  • 


(1.36) 


If  we  have 


I  =  0  » 


(1.37) 


as  would  be  the  case,  e.  g.  g  e  L4[0,L]  ,  then  the  conditions  for 

convergence  of  (1.  36)  axe  ore  stringent  than  just  the  square  summability  of 


the  av  ,  given  by  (1.  32).  We  need 

l  <l£l2+l^- l2><  - . 


<’  ••  <*  <*  -■  <*v.*  ^ 


t:wa7 


As  a  consequence  we  can  steer  (1. 21)  from  any  finite  energy  initial  state  to  a 
dense  (in  the  energy  norm)  subspace  of  final  states,  or  vice  versa,  but  we 
cannot  steer  between  arbitrary  finite  energy  states  during  [0,  ZTj]  if  (1. 37)  is 
true.  The  case  of  boundary  control  (1.  6),  already  treated  by  the  method  of 
characteristics,  will  be  discussed  more  extensively  below.  In  that  case  the 
coefficients  in  (1.  26)  are  bounded  and  bounded  below.  The  result  is, 

in  that  situation,  that  we  obtain  the  same  result  this  way  as  by  the  method  of 
characteristics  -  given  finite  energy  initial  and  terminal  states,  there  is  a 
unique  control  u  e  L2[0,  2Tj]  steering  the  one  to  the  other. 


Case  T  >  2Ti  . 


The  main  difference  between  this  case  and  the  preced- 


iuks 


LwfcS 

2, 


k  =  1,  2,  3, 


form  a 


ing  is  that  here  the  functions  e 

Rlesz  basis  for  a  proper  subspace,  E,  of  L^[0,T]  .  The  biorthogonal 
functions  p^,  p_k  exist,  but  are  unique  only  if  we  require  that  they  lie  in 
E  —  or  we  impose  some  comparable  condition.  If  we  agree  that  p^,  p  k 
belong  to  E ,  then  any  elements 


"k*  »k+  V 


p-k  =  p-k  +  q-k 


with  q. ,  q  .  «  EiC  LZ[ 0, T ]  still  form  a  biorthogonal  set  relative  to  the 

iwVs  -i«£s 

e  ,  e  .  The  convergence  properties  of  series  involving  the  pfc , 

p_k  are  much  the  same  as  in  the  preceding  case.  As  a  result  we  have  the 
same  control  capability  as  in  the  case  T  =  2T,  but  controls  are  not  unique  . 
Indeed,  if  u  is  a  control  steering  between  two  given  states,  the  family  of 
controls  u  +  r»,  u  e  E*,  all  realize  the  same  control  objective.  Again, 
this  non- uniqueness  should  be  compared  with  the  similar  property  observed  for 
T  >  2Tj  in  applying  the  method  of  characteristics. 

Using  the  theory  of  distributions  and  related  material,  boundary  value 
control  situations  such  as  (1.  6)  can  be  included  in  the  same  framework  as 
(1. 26)  but  with  g  in  a  larger  space  than  L  [0,  L]  ;  g  should  be  a  linear 
functional  (in  general  unbounded  on  LZ[0,  L]  )  whose  domain  includes  the 
domain,  A(L),  of  the  self  adjoint  operator  L»  given  by  (1.23),  with  the 
given  homogeneous  boundary  conditions.  The  gk  are  the  values  which  g 
assumes  at  the  eigenfunctions  ?k  e  6  (L)  .  A  detailed  study  of  these 
"admissible  input  elements"  is  provided  in  [17]  .  In  this  way  a  unification 
of  the  boundary  and  distributed  control  cases  may  be  achieved.  One  consequence 
of  this  is  that  the  biorthogonal  functions  pk,  p_k  which  play  such  an 


important  role  in  the  method  based  on  the  moment  problem  (1.  30),  (1. 31)  can 
actually  be  obtained  through  the  more  constructive  method  of  characteristics  as 
controls  steering  from  a  zero  initial  state  (say)  to  final  states  constructed 
using  a  single  eigenfunction  L  . 

We  began  our  discussion  here  with  the  Euler  beam  equation  (1. 1)  .  For 
definiteness,  let  us  add  a  distributed  control  term  (scalar  input)  and  specific 
boundary  conditions  so  that  we  have 


P(X) 

a2w 

+  4(EI(X,  A- 

)  =  g  (x)  u  (t) 

(1.  38) 

at2 

ax*- 

9x 

82w 

(o.t) 

83w 

(0,t)  = 

0. 

(1.  39) 

ax2 

8x3 

82w 

2 

(L,t) 

a3w 

"  3 

(L,t)  * 

0  . 

(1.  40) 

ax^ 

dx6 

In  1969  I.  P.  Quinn,  in  his  doctoral  thesis  [  27]  ,  studied  the  controllability 
properties  of  a  class  of  systems  including  this  one.  Here  the  operator 


Aw=  (EI(k)-2%.) 

P  w  8X*  8x 

on  the  domain  in  H4[  0,  L]  consisting  of  functions  obeying  boundary 
conditions  conformable  with  (1. 39),  (1. 40)  has  eigenfunctions  ^(x)  forming 
an  orthonormal  basis  for  L2  [  0,  L  ]  and  the  corresponding  eigenvalues  X^ 

grow  like  k4  as  k  —  «  ?  With  w  ^  = 

to  (128),  using  a  transformation  like  (I.  27)  applied  to  the  second  order 
differential  equations  resulting  from  the  eigenfunction  decomposition: 


X£  we  obtain  a  system  similar 


^k  +Xkwk  *  gku k  ■  0,  1,  2,  3,  ... 

( a  slight  modification  of  (1.27),  (1.28)  is  necessary  for  XQ  »  0  ;  see  [34]). 
Again  there  results 


and  the  energy  expression  for  (1.  38)  is,  equivalently. 


i  fL(p(*>  (ir-)2  +Ei(st){-^f  )2]dx 
0  az  tar 


or 


i  E  C(*fc)2  +  \(wk)2] 

k=0 


i  l  i\r  +  ickr. 


all  <  •  for  "finite  energy"  states.  ,  _  4 

lUj.  s  -lu.  s 

Quinn  was  able  to  show  in  this  case  that  the  functions  e  ,  e 

lUkS  -iyuS 

k  =  0,1,2, 3,  •••  (for  k  =  0  replace  e  ,  e  by  1,  s)  are 

2 

linearly  independent  In  L  [ 0, T  ]  for  every  T  >  0  (this  result  by  itself 


had  already  been  obtained  much  earlier  by  Ingham  [18]  who  shows ,  in  effect, 

..1 .  A.1 _ X Li X _ _  n, _ 1 a  -  x _  I 1  _ _ _ _  X  »2r«  ml  X 


that  these  functions  form  a  Riesz  basis  for  a  closed  subspace  of  L  [o,T]  for 
every  T  >  0 )  and,  additionally,  that  there  is  a  positive  number,  M(T  ), 
such  that  if  the  (non-unique)  biorthogonal  functions  p^(s),  p  ,  (s)  are 
appropriately  selected  in  L2[0,T]  ,  these  functions  are  continuous  and 
satisfy  the  point  wise  bounds 


The  fact  that  the 


-iuj'.S 


(1.  42) 


I  Puts)  |  *  M(T ),  |p  .  (s)|  *  M(T),  s  e  [0,T]  .  (1.42) 

luks  -iti»kS 

ire  e  ,  e  form  a  Riesz  basis  for  a  closed  subspace 


of  L  [0,T]  ,  T  >  0,  implies  that  initial  states  and  terminal  states  with 


(in  terms  of  (1. 41))  expansion  coefficients  0  ,  Q  and  j,  j 

can  be  steered,  one  to  the  other,  during  [0,T]  with  ue  LZ[0,T  j  , 


Vi  ,2 


^k.  1  ,2, 


provided  that 

_  n.  2  n  2  Hu  i  2  (v  <  2 

En-M-i  +i-%2-i  )<-.  E  d-iH-i 

k=0  *  *  k=0  ^  ^ 

The  boundedness  property  (1. 42)  shows  we  can  also  control  states  for  which 


I  +|3U-|  ><„,  l  (l-SU-l  +1^-1  X 


g  «  '  *  u  g. 

K  k=  0  * 


E  a- 

k=0 


^-1  +h%2-D  < 


k=o  k  * 

this  being  possible  with  a  control  function  u  (t)  uniformly  bounded  and 
continuous  on  [0,T]  . 

We  have  noted  in  connection  with  the  Timoshenko  beam  system  (1.  2),  (1. 13) 


(1. 17),  (1. 18),  that  an  adequate  control  theory,  based  on  the  method  of  character¬ 
istics,  exists  when  we  have  two  separate  control  functions,  u^(t)  and  u2(t), 
with  which  to  control  the  lateral  deflection  and  shear  deformation  separately. 

An  open  question  is  the  adequacy  of  control  using  a  single  control  input,  so 
that  (1.17),  (1.18)  becomes,  e.  g., 

jjr  (L, t)  =  au (t),  4* ( I* t )  -  (L,t)  =  &u(t) 

z  z 

with  <1  +  P £  0  .  This  problem  is  a  special  case  of  the  more  general 
question  of  the  controllability  of  linear  hyperbolic  systems  of  dimension  n  =  2m, 
involving  m  pairs  of  characteristics,  each  pair  describing  a  given  wave  mode 
propagating  in  two  opposite  directions,  by  means  of  fewer  than  m  control 
inputs.  Some  work  has  been  done  in  this  direction  by  R.  G.  Teglas  in  his 
thesis  [45]  and  by  N.  Wick  [47]  ,  but  it  is  safe  to  say  that  no  very  general 
criteria  for  this  problem  have  yet  appeared.  Particularly  valuable,  it  seems 
to  this  author,  would  be  a  study  of  the  Timoshenko  beam  system  from  the 
singular  perturbation  standpoint,  elucidating  the  behavior  of  solutions  and 
controllability  properties  as  the  modulus  of  elasticity  in  shear,  k  (x)  in  (1.  2), 
(1.  3),  tends  to  infinity. 

STABILIZATION,  CANONICAL  FORMS,  EIGENVALUE 
PLACEMENT,  etc. 

As  all  practicing  engineers  will  know,  controllability  in  itself  is  rarely 
the  prime  goal  of  control  system  design.  Stability,  and  related  criteria  such  as 
robustness.  Insensitivity  to  particular  input  frequency  bands,  etc. ,  are  more 
commonly  uppermost  in  mind.  Additionally,  there  is  the  question  of  state 
estimation  from  lower  dimensional,  noisy  observations  in  order  to  implement 
linear  feedback  control  policies.  These  subjects  have  been  pursued  almost  ad 
nauseum  for  linear,  finite  dimensional  systems.  In  the  case  of  distributed 
parameter  systems,  and  hyperbolic  systems  in  particular,  the  literature  on  this 
subject  remains  rather  sparse  and  spotty  in  its  coverage. 

As  in  the  case  of  linear  finite  dimensional  systems,  stability  and  stabili¬ 
zation  studies  for  linear  partial  differential  equations  have  tended  to  cluster 
around  two  dominant  approaches:  the  Liapounov  approach,  primarily  carried  out 
in  connection  with  systems  involving  some  form  of  "conservation  of  energy"  law, 
and  the  spectral  approach,  determining  if,  or  making  certain  that,  the 


eigenvalues  of  the  system  lie  in  an  appropriate  subset  of  the  left  half  plane.  The 
spectral  approach  suffers  from  the  disadvantages  of  greater  intricacy  of 
computation  and  the  need  to  show  that  the  spectrum  location  does,  in  fact, 
determine  the  asymptotic  behavior  of  the  system.  The  latter  brings  in  questions 
of  completeness  and  linear  independence  of  the  eigenvectors  of  the  system. 

We  will  begin  with  a  short  discussion  of  what  has  been  done  with 
Liapounov  methods.  On  the  theoretical  side  one  can  start  with  a  system 

x  =  Cx,  (2.1) 

C  generating  a  strongly  continuous  semigroup  S(t)  in  the  Hilbert  space  X 
( we  may  have  started  with  a  control  system  x  *  Ax  +  Bu,  set  u  =  Kx,  then 
C  s  A+BK)  .  We  set  up  a  quadratic  functional 

V(X)  a  (X,  QX)  , 

where  Q  is  a  bounded,  positive,  self  adjoint  operator  on  X  with  Q  >ql 
for  some  q  >  0 ,  to  serve  as  a  liapounov  function.  One  may  then  show  that 

f«r  t2  >  tj  and  x(t)  =  S  (t)xQ  a  "solution"  of  (2.1),  that 

(x(t2),  Qx(t2))  -  (x(tx),  Qx(tj)) 

f  2 

a-/  (x(s),  Wx(s))ds 

*1 

for  some  positive  self  adjoint  operator  W  so  that,  in  some  sense  which  one 
needs  to  make  precise  in  individual  cases, 

C*Q  +  QC  +  W  *  o  ,  (2.2) 

the  Liapounov  operator  equation,  is  satisfied.  An  important  result,  due  to 
Datko  [  10  ]  ,  states  that  if 

/  (x(t),  Qx(t))dt  <  -  (2.3) 

0 

for  every  initial  state  xQ  e  X,  then  the  semigroup  S(t)  is  exponentially 
damped,  i.  e. 

||S(t)H  s  Me"Yt  ,  t  a  0  , 

for  positive  numbers  M,  y  .  The  condition  (2.3)  is  satisfied  if  W  a  wl 
for  some  w  >  0  ,  as  may  easily  be  verified. 

Consider  the  linear  symmetric  hyperbolic  system  in  L*[0,  L]  = 

(L2[0,L]  )“ 

E(x)|2  =  A(x>^£  +  B(x)w+t(x,t)  (2.4) 

where  E(x),  A(x),  B(x)  are  continuously  differentiable  mXm  matrices 
defined  for  x  e  [  0,  L  ]  ,  E  (x)  symmetric  and  positive  definite,  A  (x) 


symmetric.  The  wave  and  Timoshenko  equations  can  be  written  in  this  form. 
The  "energy"  usually  is  expressed  as 

i  rL 

&(t)  =  if  (w(x,t),  E  (x) w  (x,  t))dx  , 

0 

With  appropriately  "conservative"  or  "dissipative"  boundary  conditions  at 

x  *  0,  x  =  L  ,  one  finds  that  for  t,  >  t. 

t2  L  d  1  ^  . 

C(t2)  -  e(tj)  ■  /  (if  (w(x,t),  [B(x)  +  B(x)^-A(x)]  w(x,t))dx 


+  /  (w(x,t),  f  (x,t))dx  }  dt  . 

0 

n*  * 

If  B  (x)  +  B  (x)  -  A  (x)  is  uniformly  negative  definite  or  if  the  n  dimensional 
control  function  f  (x,  t )  may  be  arbitrarily  specified  as  a  function  of  x  and 
t.  one  may  use  feedback 

f  (x.  t )  =  K  (x)  w  (x,  t )  (2. 5) 

in  such  a  way  that  . 


2  L 


C(t2)  -  eftj)  =  -f  f  (w(x,  t),  W  (x)  w  ( x,  t )) dx  dt 


(2.  6) 


with  W(x)  uniformly  positive  definite  and  symmetric  on  [0,  L]  .  Then  one 
can  apply  Datko's  result,  or  more  simple  arguments,  to  show  that  solutions  of 
(2.3),  (2.4)  are  uniformly  exponentially  damped  in  I^[0,  L]  norm. 

Note,  however,  that  if  B  (x)  +  B  (x)*  -  A  (x)  =  0  or  for  some  other 
reason  falls  to  be  positive  definite,  and  If 

f(x,t)  =  D  (x)u  (x,  t)  or  f(x,t)  =  D  (x)  u  (t) 
with  dim  u(x,  t)*r<R  in  the  first  instance,  u  a  function  of  t  only  In 
the  second  Instance,  then  we  cannot,  in  general,  achieve  (2.5)  with  W(x) 
uniformly  positive  definite.  Comparable  difficulties  arise  when  boundary 
control  is  employed.  In  such  cases  it  is  a  form  of  the  La  Salle  "invariance 
principle"  (see,  e.  g.  [19])  which  must  be  appealed  to,  rather  than  the  basic 
Liapounov  theory,  for  an  analysis  of  presumed  asymptotic  stability  properties  of 
the  system.  This  has  been  discussed  in  some  detail  in  [36]  and  [33]  and 
we  give  only  the  briefest  outline  here. 

The  "invariance  principle",  as  it  applies  to  finite  dimensional  systems, 
relies  heavily  on  the  compactness  of  the  "  w-limit  set"  of  the  system  in  order  to 
reach  the  final  conclusion  of  asymptotic  stability.  Comparable  compactness 
properties  associated  with  the  solutions  of  an  infinite  dimensional  system  are 


generally  difficult  to  realize  but  the  initial  attempts  to  extend  the  theory  nonethe¬ 
less  relied  on  establishing  some  sort  of  compactness  property.  One  of  the  first 
contributions  in  this  direction  was  due  to  Dafermos  [  9  ]  who  studied  weak 
damping  of  the  wave  equation,  relying  on  the  almost  periodic  nature  of  the 
system  solutions  to  provide  the  required  compactness.  Slemrod  [  43  ]  studied 
the  boundary  damped  wave  equation  by  introducing  suitably  weakened  topologies 
—  as  compared  with  the  usual  topology  associated  with  the  energy  norm  —  and 
was  able  to  conclude  a  correspondingly  weakened  form  of  asymptotic  stability. 

Knowing  that  controllability  implies  stabilizability  in  the  case  of  autono¬ 
mous  finite  dimensional  linear  systems,  we  are  not  surprised  to  find  control¬ 
lability  playing  a  role  in  the  study  of  asymptotic  stability  and  stabilization 
properties  of  autonomous  infinite  dimensional  linear  systems.  This  is  discussed 
in  some  detail  in  the  paper  [28]  by  J.  P.  Quinn  and  the  author  and  also  in 
[33],  Systems  of  the  form  (2.3),  but  with  the  control  appearing  in  the  bound¬ 
ary  conditions,  are  studied  in  [28]  prior  to  the  main  discussion  on  the  boundary 
damped,  higher  dimensional  wave  equation.  We  can  give  an  idea  of  the  flavor 
of  the  arguments  employed  using  a  simple  example  based  on  the  wave  equation 
(1.  4)  with  control  appearing  in  the  boundary  conditions  (1. 5),  (1.  6).  If  in  this 
system  one  employs  the  feedback  law 

«(»)  *  -Y  X  (2.7) 

the  closed  loop  system  is  (1.4),  (1.5)  together  with  the  "closed -loop"  Robin 

type  boundary  condition 

<i.t)  +  YTjJr(l.t)  *  0  .  (2.8) 

Here  a  short  computation  shows  that  with  the  energy  £(t )  defined  by  the 
expression  (1.7)  we  have,  for  t2  >  tj, 

*2 

e(t2)  -  eftp  =  -yp  (L)  /  ^-(L,t)2dt  .  (2.9) 

It  Is  not  feasible  to  fit  this  situation  into  the  general  pattern  based  on  the 
Llapounov  operator  equation  (2.  2)  but,  since  we  expect  (correctly)  that,  along 
with  (1.4),  (1.5),  (2.8) 

(L,t)  =  0  =$>  w(x,t)  =  0, 

an  "invariance  principle"  type  of  argument  appears  to  be  in  order.  But  we 
will  use  a  variation  on  this  procedure  which  makes  use  of  the  controllability 
already  established  in  Section  1,  Let  v(x,t)  be  a  controlled  solution  of 


7*7*5  V-  ’>  \*  ^  r  *■..»  ..■  .-.  V-  ^  , ! 


(1. 4),  (1.  5),  (1.  6),  u  (t)  being  selected  so  as  to  steer  the  initial  state 

v(x,0)  a  w(x,0),  |jr  (x,0)  =  ■—  (x,0),  (2.10) 

agreeing  with  the  initial  state  of  the  solution  w  (x,  t)  of  (1. 4),  (1. 5),  (2.  8), 
to  the  zero  final  state 

vfx,^)  =  0,  —  (x.Tj)  =  0.  (2.11) 

Tj,  as  described  earlier.  Defining  the  "energy  inner  product" 

L 


(w(-,  t),v{*  ,  t))  =  f  [p  (x) (x,t)  (x,t)  +  p(x)  |^(x,t)|~(x,t)]  dx 

it  is  found,  using  (1. 5),  (2.  8),  (2. 10),  (2. 11),  that 

<w(.,0),  v(-^0)>  -  (w(* ,  2Tj),  v(«,  2Tj)>  =  ||w(.,  0)|||  = 

s  “P<L>^  X[|^(^t)|f(L,t)+g(L,t)^-(L,t)]dt 

*  pjl)/^1  ~r  [v|f  (L,t)  +  ~  (L,t)]  dt 
2T 

«p<uf  1T?(Ut)tv  +u(t)]dt. 

Here  H  w ( ■  ,  0  )[| ^ ,  the  energy  norm  at  t  =  0,  is  2fi(0)  .  Applying  the 
Schwartz  inequality 


4C(0)' 


.ZT 


p(l)jT"1  (L,trdt 


.2T. 


(2.12) 


fQ  1(Y^(I',t)+u(t))2dt. 


A  slightly  more  detailed  study  of  the  control  problem  for  (1. 4),  (1.  5),  (1.  6)  in 
the  case  T  s  21,  (or  T  >  2T, )  shows  that  control  from  an  initial  state 
w(x,  0),  (x,  0)  to  0,0  at  time  2T^  is  realized  with  a  control  u(t) 

which  satisfies  2? 

/  u(t)ZdtsLe(0) 

0  0 

and,  for  the  resulting  controlled  solution  we  have 

j“‘  |2(L,t)2dt  s  K^iO) 

for  certain  positive  constants  KQ  and  .  Then  (2. 12)  easily  yields 

YP (L)f  (L,t)2dta  - - 

0  31  2(KQ+YZK1)e(0) 


2  Ke(0) 


and,  setting  tj  =  0  ,  t2  =  ZI^  in  (2.9),  we  have 

e(2T1)  S  e (o)  -  Ke(o)  =  (l-K)e(o).  (2.13) 

Since  g(2T^)  is,  from  (2.9),  (2.13),  positive  and  less  than  or  equal  to  £(0) 

we  conclude  0  <  1  -  K  <  1  . 

Repeating  the  above  argument  on  successive  intervals  [0,  2T^]  ,  [2T^, 

4TX]  •••  [2kTj,  20^+1)^]  , . . .  and  using  the  monotonicity  of  £(t),  as 
implied  by  (2.9),  we  conclude  that  e(t)  decays  exponentially  to  0  as  t  — «. 
The  same  general  argument  can  be  used  with  a  fairly  wide  class  of  boundary 
damped  linear  symmetric  hyperbolic  systems  (2.  4)  and  with  many  other  systems 
which  are  energy  conserving  in  the  uncontrolled  situation  and  suitably  strong 
controllability  properties.  The  Timoshenko  system  (1.2),  (1. 3),  with  appropriate 
boundary  conditions,  is  in  this  class.  As  far  as  the  author  is  aware,  the  Euler 
beam  model  (1. 1)  has  not  yet  been  studied  from  this  point  of  view. 

The  spectral  approach,  as  we  have  already  indicated,  involves  a  direct 
analysis  of  the  eigenvalues  and  eigenfunctions  or,  more  generally,  the  spectrum 
and  invariant  subspaces,  of  the  generating  operator  C  for  a  given  system 
x  =  Cx,  possibly  derived  from  a  control  system  x  =  Ax  +  Bu  by  the  use  of 
linear  feedback  u  =  Kx  so  that  C  =  A  +  BK  .  A  fairly  common  case,  which 
can  be  treated  with  minimal  difficulty,  arises  when  all  but  finitely  many  of  the 
eigenvalues  of  C  have  negative  real  parts.  Under  genetically  valid  control¬ 
lability-type  conditions  it  is  then  possible  to  move  the  unstable  eigenvalues 
into  the  left  half  plane  while  either  keeping  the  stable  eigenvalues  fixed  or  else 
maintaining  a  certain  margin  of  stability.  Work  of  this  sort  has  been  carried  out 
byTrlggiani  [46]  ,  Sakawa  [40]  ,  [41]  and  others. 

A  somewhat  more  challenging  task  arises  when  one  starts  with  a  system 
having  infinitely  many  eigenvalues  in  the  closed  right  half  plane  (usually  one 
considers  a  conservative  system  wherein  all  of  the  eigenvalues  of  C  are  purely 
imaginary)  and  one  attempts  to  devise  a  feedback  law  to  move  all  of  these 
eigenvalues  over  into  the  open  left  half  plane.  A  number  of  procedures  have  been 
examined  in  this  connection. 

In  [32]  a  second  order  system  with  scalar  control 

x  +  Ax  =  bu,  x,  b  e  X  ,  (2. 14) 

is  studied,  X  being  a  real  Hilbert  space  and  A  an  unbounded  positive  self 
adjoint  operator  on  X  .  Assuming  that  A  has  a  Riesz  basis  of  eigenvectors 
k=  1,  2, 3,**«,  in  X,  and  corresponding  positive  eigenvalues  , 


(2.15) 
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increasing  with  k  ,  ks  1,2,3, 


x  and  b  may  be  expanded  as 


x=  E  Vk>  b  =  E  Mt 


k=l  k=  1 

convergent  in  X  ,  with  square  summable  coefficients, 
minimal  condition  for  approximate  controllability 


We  assume  the 


K  t  0  ,  k  =  1,  2,  3, 


•  • 

The  energy  form  is  £  [  (x,x )  +  (x.  Ax )]  =  g  and  elementary  computations 


show  that  for  (2. 14)  and  for  any  T  >  0 


e(T)  -  £(0)  =  /  (x(t),b)u(t)dt  . 
0 


(2.16) 


(2.17) 


It  follows  that  with 

u(t)  =  -y(x(t),b) 
the  energy  C(t)  is  non- increasing  with  increasing  t.  So  far  this  is 
basically  a  Liapounov  approach  employing  what  is  known  in  the  engineering 
literature  as  an  ILAF  ( Identical  Location  of  Accelerometer  and  Forces) 
approach.  The  resulting  closed  loop  system  is,  still  in  second  order  form, 

x  +  Bx  +  Ax  =  0  (2. 18) 

with  B  defined  by 

Bx  =  y  (x, b)b  .  (2.19) 

With  y  =  x ,  one  may  consider  the  equivalent  first  order  system  in  XX  X, 


(*)-  ■(:)■.•■£  ■•)■ 


(2.  20) 

and  ask:  what  are  the  eigenvalues  and  eigenvectors  of  C?  It  is  here  that 
one  leaves  the  second  method  of  Liapounov  and  returns  to  his  first.  In  [32] 
a  perturbation  analysis  is  carried  out,  valid  for  small  values  of  y  in  (2. 17), 
(2.19).  It  is  shown  that,  under  the  separation  assumption 

wk+1  -  wk  >  d  >  0,  ,  (2.  21) 

the  eigenvalues  of  C,  which  for  y  =  0  are  ±iw^ ,  k  -  1,  2, 3,  •  •  • ,  all 
have  negative  real  parts  for  y  >  0  and,  moreover,  designating  the  perturbed 


eigenvalues  by  ^( y), 
we  have  (cf.  (2. 21)) 


k  =  ±  1,  +  2,  ±3, 


Ck(0>  =  iwk'  C-k(0>  =  "iwk 


Ck(Y)  =  iwk  -  7  \\\2  +  ®(Y2' 


).  k  - 


(2.  22) 


C_k(V  )  =  -iw,  -  j  |  b,  p  +  ®  (y 


l  wvT 


), 


It  is  also  possible  to  show  that  the  perturbed  eigenvectors  cont 
Riesz  basis  for  the  space  XXX.  From  this  it  follows  that 
(2. 18)  tend  strongly  to  zero  in  the  energy  norm,  though  not  at  < 
ential  rate. 

Following  Wonham' s  initial  results  [48]  on  the  finite  c 
there  has  been  considerable  interest  displayed  in  the  question 
determination  via  linear  feedback  for  distributed  parameter  sysl 
of  the  system  (2.14),  equivalently, 

O  I 

-A  O 

with  initial  (u  =  0)  eigenvalues  ±  iw.'  k=  1,  2,  3,  ••• 
be  phrased  as  follows:  we  suppose  use  of  a  linear  feedback  f 

u  =  (A2x ,  k^)  +  (y,k2),  kj,  k2  e  X, 

i 

bounded  relative  to  the  energy  norm  (x,Ax)  +  (y,  y)  =  (A2x, 
XXX.  With 

KjX  =  (A2x,  kx)b,  K2y=  (y,k2)b 
the  closed  loop  system  is 


u. 


One  can  now  ask:  What  eigenvalues  can  be  achieved  for  the  < 
system  (2.24)  by  appropriate  selection  of  k^,  k2  in  (2.  2E 
time  the  author  was  under  the  impression  that  his  approach  via 
[35]  (more  on  this  below)  was  the  first  treatment  of  this  qi 
fact,  it  appears  that  this  credit  must  go  to  Prof.  Sun  S.  -H.  of 
University  who  treated  this  problem  by  a  more  sophisticated  a* 
perturbation  technique  used  by  the  author  in  [32]  to  obtain  t 
Sun  was  able  to  show,  with  an  assumption  similar  to  (2.  21)  a 
basis  assumption  on  the  open  loop  eigenvectors,  that  the  total 
achievable  by  use  of  (2.  25)  coincides  with  sequences 
1,  2,  3,  •  •  •  for  which,  assuming  the  b.  ^  0  as  before. 
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His  very  important  paper  has  been  translated  by  Ho  L.  -F.  in  [44]  .  Some 
comparable,  but  necessarily  weaker,  results  have  been  obtained  by  Reid  in  his 
thesis  [29]  for  the  equation  of  linear  surface  waves  where  (2.21)  is  not 
satisfied  and,  in  fact,  lim  (w^+^  -  =  0  .  Other  results  in  this 

direction,  for  hyperbolic  systems  of  various  types,  have  been  obtained  by  Clark 
[5],  [6]  and  by  Ho  in  his  thesis  [16]. 

Much  of  the  initial  impetus  for  the  study  of  control  canonical  forms, 
both  for  finite  and  infinite  dimensional  systems,  came  from  the  spectral 
determination  question  discussed  above,  but  the  subject  is  interesting  in  its 
own  right  and  shows  some  promise  of  being  adaptable  for  "real  world"  control 
implementation.  The  reader  will  recall  that  a  finite  dimensional  controllable 
system 

x  =s  Ax  +  bu,  x  e  Rn , 

with  scalar  control  u  is  equivalent,  via  a  state  space  similarity  transform¬ 
ation  (see  [20],  [35])  to  a  system  in  rational  canonical  form  corresponding 
to  the  n-th  order  scalar  equation 


y(">  +  a.y'"-1’  + 


where 


y  +any  =  u. 


(2.  27) ' 


p(X)  =  det(XI-A)  =  \n  +  a,  Xn_1  +  •••  a  A  +  a 

•  n**i  n  * 

is  the  characteristic  polynomial  of  the  matrix  A.  Comparable,  but  somewhat ! 
more  intricate,  results  are  available  for  systems  with  higher  control  dimension  ! 
[20]  ,  [2]  .  In  [38]  we  note  that  if  one  employs  a  scalar  linear  observation  j 


y  =  h  x  =  ( x,  h ), 


(2.  28) 


there  is  exactly  one  observation  vector  he  Rn  for  which  (2.  28)  satisfies  ! 
(2.  27);  for  general  h  the  right  hand  side  will  involve  the  derivatives  of  u  1 
of  order  s  u  -  1  .  Systems  (2.  27)  are  particularly  easy  to  deal  with. 

Closed  loop  eigenvalues  ^ *  **,  £n  may  be  realized  simply  by  forming 
the  polynomial 


q(M  =  i  I  -  Cic)  =  xn+ci 

k  =  1  *  1 


\n~l  + 


+  c_  ,  x.  +  c 
n-i  n 


and  determining  u  by  linear  feedback  on  the  observation 


and  its 
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derivatives , 


u  ■  I  <ak-ck>y 


(n-k) 


Apparently  less  well  known,  but  quite  obvious,  is  that  the  control  problem  for  i 

j 

(2.  27)  is,  in  a  sense,  trivial.  Let  us  suppose  the  initial  instant  is  taken  to  i 
be  t  s  0  and  control  is  to  be  effected  during  0  s  t  s  T  .  Let  the  initial! 
state  be  specified  by 


y<n"k)  (0)  =  Vk+1, 


and  the  terminal  state  by 


k  =  1,  2, 


y'“"~'  CT)  =  Yn  k+1  ,  k  =  1,  2,  •  •  • .  n  .  (2.  30)  j 

If  y(t)  satisfies  (2.29)  and 

y<n>(t)  =  v(t),  0  stsT  (2.31) 

then  we  see  readily  that  for  k  =  1,  2,  •  •  • ,  n 

y(n-k>(t).^  yn_1+1  +  ^(t-sf-Mslds  . 

and  (2.  30)  is  achieved  just  in  case 

f  <t-s)k_(  v(s)ds  .  Yn_k.1  -  £  yn-J-l  TkTTp.  • 

U  £  =  1 

k  =  1,  2,  • .  • ,  n  . 

This  is  easily  solved  for  v  in  various  function  classes,  e.g.  polynomials  of 
degree  ^  n  -  1,  etc.  and,  it  should  be  noted,  the  solution  has  nothing  to  do  I 
with  the  coefficients  in  (2.  27)  so  the  calculation  can  be  carried  out  once  for 
any  given  T  and  recorded  for  use  ever  after.  Then  in  a  given  canonical 

i 

system  (2. 27)  we  need  only  set 

n  j 

u(t)  =  v(t)  -  Yj  akyln_K,(t)  (2.  32)  i 

k=l  | 

to  realize  the  desired  control  objective. 

Since,  in  a  given  control  context,  it  is  not  likely  that  the  available 

observation  (2.  28)  will  be  the  particular  one  for  which  (2.  27)  obtains,  the 

above  result  might  seems  to  be  a  generally  useless  curiosity.  It  turns  out, 

however,  that  in  canonical  form  theory  there  is  a  counterpart  to  the  more  widely 

known  observer  theory.  If  C  is  any  nXn  matrix  whose  minimal  and 


(2.  29) 


(2.30) 


(2.  31) 


b  "j.  v  v.  *>  »>  ■>  n-v  v«i  »jt  «i  »>  »;  ■;»,«»  p « jnj.  v'1-  ,.i,.»,.i «.« «.i  ■>  -.i 


characteristic  polynomials  coincide,  it  is  possible  to  select  (non-uniquely)  r, 
d  and  J  such  that  the  augmented  system 
x  =  Ax  4-  bu 


z  =  ry  +  Cz  +  du  (=  rh  x  +  Cz  +  du  since  y  =  h  x) 

with  augmented  observation 

w  =  y  +  Jz  =  h"x  +  J^z 

is  in  canonical  form,  so  that  for  some  coefficients  o^,  a 2,  •  •  • , 
w(2n)  +  a^w(2n-l)  +  . . .  +  <*2m.  1  w'  +  a  2mW  =  U  * 


(2.  33)  j 
(2.  34)  . 

(2.  35) 


2m 


The  adjoined  system  (2.  34)  can  be  realized  electronically,  just  as  an  observer: 
system  is,  and  the  considerable  freedom  in  choice  of  C  ,  r,  d  and  j 
provides  much  design  flexibility.  In  some  cases  the  dimension  of  (2.  34)  can  : 
be  reduced.  The  proof  that  (2.  33),  (2.  34),  (2.  35)  can  be  made  a  canonical  I 

system  appears  in  [38]  . 

A  parallel  control  canonical  form  theory  has  been  developed  for  certain 
hyperbolic  distributed  parameter  systems,  involving  neutral  functional  equations 
in  place  of  the  n-th  order  scalar  equation  (2.  27).  The  theory  is  quite  complex, 
especially  as  it  applies  to  partial  differential  equations  with  variable 
coefficients  (see  [35],  [16]  ,  [  38]  ,  [  39  ]  e.g. ).  To  give  an  idea  how  the 
theory  is  developed  we  will  consider  the  constant  coefficient  case  of  (1.  4) 
which,  without  loss  of  generality,  we  can  take  to  be 


32w 


at 


a2w 


ax 


=  o,  t>o,  osxsi, 


3w 


(2.  36) 


w(o.t)  =0,  ~  (l,t)  =  U  (t)  . 


(2.37) 

The  normalized  eigenfunctions  of  the  corresponding  homogeneous  system  are 

2k-l 


trx,  k  =  1,  2,  3, 


Setting 


wk  = 


^k(x)  =  '■PT  sin  — 2 
2k -1 

— 2~  *  and  forming  the  expansions 

m 

w(x,t)  =  Yj  wk(t)  ^k*x>  ' 
k=  1 

m 

lT(x't>=  E  vk(t)5*k(x)’ 

k=  1 

followed  by  the  transformation 


(2.  38) 

(2.39) 

i 

(2.40)  , 


(2.4i)  ; 


we  have,  for  k  =  1,  2,  3, 


V  =  iwu  T,.  + 


*  +  u(t)*  “  -1WA +  ^ u 


Consider  now  the  neutral  delay  equation 


(2. 42) 


y(t+2)  +  y(t)  *  u  (t  +  2 )  ,  (2.43) 

The  characteristic  function  of  the  homogeneous  equation  is 

p(M  =  e^  +1  =  2e  ^  cosh  X. 

and  the  zeros  of  p(M  are  precisely  the  eigenvalues  ±iw^  appearing  in 
(2. 42)  .  The  transfer  function  for  (2. 43)  is 

.  2 1  »  i 


T„(M  = 


e  + 1 


which  can  be  rewritten  as 


Vx>-  Z  T TT-T-  +  “k“  “T4-  " 

k=l  *  +  wk 

If  we  define  an  observation  y  (t)  on  (2. 42)  by 

Y(t>=  E  thknk(t,+gk^k(t)1  +  iu<t> 
k  =  1 

the  transfer  function  for  y  is,  formally. 


_  1  sinh  X.  ,  1_ 
If  cosh  k  2 


_  2k-l 


(2.44) 


(2.  44) 


y  hk(-I)k~X  V2 

*]_  *•  ’  luk 


7t 

^  +  iui. 


which  may  be  seen  to  agree  with  (2.  44)  just  in  case 

„  k  - « k  -  . 

*  V2 

Using  (2.38),  (2.39),  (2.41),  (2.42)  it  may  be  seen  that  this  choice  of  h^,  gj 

corresponds  to 


■  \  +  Tu«l>  -  *  <T(r(>.t>  +  15-0.0) 


(2.45) 


This  observation  on  (2.  42),  and  no  other,  satisfies  the  scalar  equation  (2. 43) 
which  serves  as  the  control  canonical  form  for  (2. 42).  The  details  of  the  above 


>  -  v>iVvV.v 


mow 


calculations  and  some  idea  of  the  form  of  a  general  theory  appear  in  (  38  ]  and 
[39]. 

If  the  canonical  observation  (2.  45)  were  actually  available,  so  that 
we  have  (2.43),  its  usefulness  is  quite  clear.  For,  with  the  causal  feedback 
law  , 


u  (t  +  2)  s  (1- v)y(t)  -  /  c (s)y  (t  +s)ds 

0 

(2. 43)  transforms  to 

2 

y  (t +2)  +  yy(t)  +  f  c(s)y(t  +  s)ds  =  0 

0  & 
and  it  is  known  from  [35],  [44]  that  the  exponential  solution  e  , 

of  (2.  47)  can  be  made  such  that 


(2. 46) 


(2. 47) 
<-k* 


Ck  =  iuk  +  a  +  V  C-k  =  iwk  T  T  &-k  * 


+  a  +  e 


where  a  is  a  complex  number  (ordinarily  negative)  determined  by  y  and 


are  arbitrary  complex  numbers,  determined  by  c  c  L^[0,  2]  ,  such 

Q» 

i 2  .  •  ,2 


ck'  -k 
that 

Z  <KI2  +  Kkh  <  -  • 

k  =  1 

It  may  be  shown  that  these  are  the  eigenvalues  of  the  closed  loop  system  (2.36), 
(2.37),  (2.45).  (2.46). 

In  a  given  application,  however,  it  is  entirely  likely  that  the  particular 
"canonical"  observation  (2.46)  will  not  be  available.  Indeed,  in  the  examp  14 
indicated,  since  this  observation  is  taken  at  the  same  point  where  control  is 
applied  and  might,  therefore,  be  subject  to  a  certain  amount  of  noise 
disturbance,  it  might  not  be  desirable  to  use  this  observation  in  practice.  To 
illustrate  the  use  of  the  technique  of  canonical  augmentation  (or  "canonical 
compensation",  perhaps)  let  us  consider  the  same  system  (2.36),  (2.37), 
but  suppose  the  available  observation  is 


y(t)  =  ^  (O.t)  . 

It  is  not  hard  to  show  in  this  case  that  y  (t)  satisfies 

y  (t  +  2)  +  y  (t)  =  u(t  +  l) 


(2.48)  j 


(2.49) 


rather  than  (2.  43).  This  "central"  control  canonical  form  is  not  as  usable  as 

the  "backward"  form  (2. 43)  because,  unlike  (2. 46), 

f2 

u(t  +  l)  *  (l-y)y(t)  -  J  c(s)  y(t  +  s)ds 

0 


is  not  a  causal  feedback  law  and  cannot  be  implemented.  But  now  couple  (2.49) 
with 


and  let 


z(t+2)  +  Pz(t)  a  au(t+2)  +  bu(t+l) 

+  cy  (t  +  1)  +  dy  (t ) 


(2.50) 


w(t)  =  y (t)  +z(t)  . 

One  ordinarily  will  take  |pj  <  1  so  that  the  homogeneous  part  of  (2.50)  is  - 
asymptotically  stable,  thus  avoiding  the  growth  of  parasitic  solutions  in  the 
compensator.  Since 

[y(t+4)  +y(t+2)  -  u(t+3)]  +  p  [y(t+2)  +y(t)  -  u(t+l)]  =0 

while 

(z  (t  +4)  +  pz  (t  +  2)  -  au  (t+4)  -  bu  (t  +  3)  -  cu(t  +1)  - cy(t  +3)  -dy  (t  +2)] 

+  [  z(t  +2)  +  pz  (t)  -  au  (t  +  2)  -  bu  (t  +1)  -  cy  (t  +1)  -  dy  (t)  ]  =  o 
we  find  that 

w(t+4)  +  (l+  p)  w  (t  +  2)  +  pw(t)  a  au  (t  +  4)  +[l+b]  u(t+3) 

i 

+  au (t  +2)  +[p  +  b]  u(t+l)  +  c[y (t  +3)  +y(t+l)] 

+  d[y(t  +  2)  +y  (t) ]  a  (using  (2.49)) 
au (t  +  4)  +  [1+b]  u  (t  +3)  +  [a  +c]  u(t  +  2)  +[p+b+d]  u(t+l)  . 

Then  it  is  easy  to  see  that  with 

a  =  l,  l  +  b=  a+  c=p+b+d  =  0, 

i.  e.  with 

a  =  1,  b  a  -1,  C  a  -1,  d  a  1  -  P  ,  ! 

we  arrive  at  the  "backward"  canonical  form  satisfied  by  w(t)  :  ! 

w(t  +  4)  +(l  +  p)  w(t+2)  +  pw(t)  a  U(t  +  4) 

for  which  causal  feedback  laws 

u  (t+4)  =  -YjW(t+3)  +  [1+p  -y2]  w  (t  +2)  -  y3  w(t  +  l) 

+  [p-Y4]  w(t)  -  /  c(s)w(t  +  s)ds  (2*51) 

0 

may  be  implemented,  yielding  overall  closed  loop  systems 


w(t  +  4)  +Y,w(t+3)  +Y,w(t  +  2)  +  Y^w  ( t  + 1)  +Y/fw(t) 

4  4  (2  52) 

+  f  c(s)w(t+s)ds  =  0. 

0 

It  is  necessary  to  check  separately  that  the  system  (2.  49),  (2.  50),  (2.  51)  is  ! 

observable  in  any  given  case. 

The  exponential  solutions  of  (2.  52),  and  hence  the  eigenvalues  of  | 
(2.36),  (2.37),  (2.48),  (2.50),  (2.51)  may  be  determined  with  the  same 
flexibility  as  already  noted  for  (2. 47).  This  is  discussed  in  some  detail  in  the 
thesis  of  R.  G.  Teglas  [45]  .  A  complete  theory  of  canonical  compensation 
for  hyperbolic  systems  remains  to  be  developed  but,  we  hope,  the  example 
given  here  gives  reason  to  believe  that  the  method  is  a  promising  one.  It  is 
clear  that  there  are  some  connections  with  observer  theory  as  developed  in 
[22]  and  elsewhere;  these  connections  remain  to  be  worked  out. 
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Admissible  Input  Elements 
for  Systems  in  Hilbert  Space 
and  a  Carleson  Measure  Criterion" 

by  E.F.  Ho+  and  D.  L.  Russell++ 

Abstract 

We  study  the  control  system 

x  *  Ax  +  bu,  x  €  X  ,  u  scalar, 

where  A  generates  a  semigroup  on  the  Hilbert  space  X  but,  in  general, 
the  control  input  element  b  X  .  Many  boundary  value  control  systems, 
point  control  force  situations,  etc. ,  can  be  studied  in  this  context  .  We 
define  and  analyze  "admissible"  input  elements  b  and  develop  sufficient 
conditions  for  b  to  be  admissible  in  terms  of  the  Carleson  measure  theorem 
of  H^-theory  . 
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1.  Introduction .  One  commonly  studies  linear,  time  invariant  control 
systems  in  a  Banach  space  X  in  the  form 

x  =  Ax  +  Bu,  x  e  X  ,  u  e  U , 

where  A  is  the  generator  of  a  strongly  continuous  semigroup  of  bounded 


(1.1) 


operators  (S(t)|  t  2=  0}  on  X  and  B  is  a  bounded  operator  from  the 
control  space,  U,  into  X.  If  u:[0,®)  — U  is  locally  (Bochner) 
integrable,  generalized  (or  "mild")  solutions  of  (1.1)  corresponding  to 


an  initial  state 


x(0)  =  xQ  e  X 


can  be  represented  by  the  "variation  of  parameters"  formula  (  see,  e.g.  [3]  , 

[11]) 

x(t)  =  S(t)x_  +  f  S  (t  -  s)  Bu(s)  ds  (1.2) 

0  0 

and  a  number  of  properties  of  x(t)  thereby  deduced  . 

It  is  well  known,  however,  that  most  of  the  "interesting"  infinite 
dimensional  control  systems  do  not  arise  this  way  because  the  degree  of 
controllability  of  a  system  (L  1)  with  B  bounded  is  rather  restricted  if,  as 
is  usually  the  case,  U  is  finite  dimensional  or  for  some  other  reason  the 
operator  B  is  compact  .  Indeed,  most  of  the  mathematically  intriguing 
examples  arise  in  the  context  of  partial  differential  equations  with  boundary 
value  control  inputs,  control  forces  exerted  at  isolated  points,  etc. ,  and  in 
the  context  of  functional  equations  which  involve  values  of  the  control  of 
discrete  instants  ,  viz. ;  u  (t) ,  u  (t  -  T^), ...  ,  u  ( t  -  Tn)  .  In  each  of  these 
cases  the  formulation  (1. 1)  is  inadequate  and  one  must  consider  input 
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operators  '  B  whose  range  is  not  restricted  to  the  space  X  . 

A  number  of  authors  have  addressed  the  problem  of  interpretation  of  (1. 1) 
for  operators  B  of  rather  general  type  .  We  particularly  cite  the  contributions 
of  Curtain  and  Pritchard  [  3  ]  ,  Zabczyk  [22]  ,  Fattorini  [  6]  ,  and  Washburn  [  20]  . 
It  seems  fair  to  say  that,  as  brought  out  in  [  3  ]  ,  the  theory  is  more  extensive 
and  generally  applicable  in  the  case  of  systems  of  "diffusion  type",  ordinarily 
involving  holomorphic  semigroups,  than  in  systems  of  "wave"  or  hyperbolic 
character  . 

In  the  present  article  we  shall  restrict  our  attention  to  spaces  X  which 

are  separable  Hilbert  spaces  and  to  finite  dimensional  control  spaces  U  . 

Taking  U  to  be  Rm  ,  (1. 1)  becomes 

m 

x  =  Ax  +  £  bj  uJ  (1. 2) 

j=  1 

where  b^  is  the  control  input  element  associated  with  the  j-th  control 
component  .  Since  every  solution  of  (1. 2)  is  a  linear  combination  of 
solutions  of  x  =  Ax  and  the  individual  systems  x  =  Ax  +  b^  u^  , 
j  =  1,  2, . . . ,  m  ,  we  may,  without  loss  of  generality,  confine  our  discussion 
to  systems 

x  =  Ax  +  bu  (1.  3) 

wherein  the  control  u  is  scalar  valued.  Much  of  our  theory  can  be 
extended  to  cases  wherein  U  is  infinite  dimensional  but  we  will  not  do 
that  here  . 

What  distinguishes  the  present  study  from  earlier  contributions  is  the 
attention  which  we  pay  not  only  to  the  relationship  between  the  operator.  A, 
and  the  input  element  b  ,  but  also  to  the  relationship  between  b  and 
the  semigroup  S(t)  generated  by  A  .  In  cases  where  A  has  discrete 


spectrum  {\^  J  k  e  K}  ,  K  being  a  countable  index  set, 
to  a  study  encompassing  the  input  element  b  ,  the  eigenve 
of  A  ,  the  corresponding  eigenvectors  of  the  dual  operator, 
defined  in  Section  2  ,  and  the  exponential  functions  exp  \ 
It  is  in  particular  reference  to  the  latter  that  what  is  probably 
important  idea  of  this  paper  is  developed  .  We  show  that  a 
condition  for  b  to  be  an  "admissible  input  element"  (del 
can  be  given  in  terms  of  a  measure  on  Borel  subsets  of  the  cor 
whose  support  is  {-kjjk  e  K}  .  When  that  measure  turns 
Carleson  measure  the  input  element  b  is  admissible  .  Thi: 
out  yet  again  the  intimate  relationship  between  the  control  the 
dimensional  linear  systems  and  parallel  developments  in  Hp 
([5],  [8],  [12]  )  and  the  related  theory  of  completeness  ai 
of  sets  of  complex  exponentials. 
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2.  Admissible  Input  Elements  . 

Let  X  be  a  separable  Hilbert  space  and  let  A  be  a  closed  operator 
on  X  with  domain,  g  (A)  ,  dense  in  X  ,  generating  a  strongly  continuous 

semigroup  of  bounded  operators  S(t)  on  X  for  t  2:  0  ,  For  b  e  X 

the  (generalized,  or  "mild")  solution  of 

x  =  Ax  +  bu ,  U€Liocf°'"^'  f2*1) 

x(0)  *  xQ  e  X  ,  (2.  2) 

is  given  by  the  "variation  of  parameters"  formula 

x(t)  =  S(t )xQ+J  S(t-  s)bu(s)ds  ( ’••  3) 

and  may  be  seen  to  be  a  continuous  function  x  :  [  0,  ®)  X  .  Whether 
x(t)  is  defined  for  each  t  a  0  and  (2.1)  holds  is  more  complicated: 
sufficient  conditions  are  that  b  €  fi(A)  or  that  u  is  differentiable  as 
a  function  of  t  ([3],  [11]). 

In  this  paper  we  wish  to  consider  (2. 1)  ,  (2.  2)  in  certain  cases  where 

b  does  not  lie  in  X  and  to  provide,  for  such  b,  a  formula  parallel  to 

(2.  3)  .  Our  approach  is  similar  to  that  used  in  [14]  . 

Identifying  X  with  its  dual  X  ,  we  denote  the  duality  relationship 
by  (x,  y)  ,  x  c  X ,  y  e  X  ,  linear  in  both  x  and  y  .  Where  X  is 

the  complexificatlon  of  a  real  Hilbert  space  X^  the  conjugate  element  *y 

is  well  defined  for  each  y  e  X  and,  with  (  ,  )  denoting  the  inner 

product  in  X , 

(x,  y)  =  <x,y>  ,  (x,y>  =  (x,y)  . 

The  bilinear  form  <  ,  )  is  symmetric,  i.  e. ,  (x,  y)  =  <y,x>,  x,y  eX, 


v'  •/  N  V* 


.  •  ‘ . 


w 


5 

r~- 

a 


and,  for  all  xeX, 


|<x,y>| 

|jxjjx=  sup  - 

x  ycx  llvllx 

y  ^  0  * 


(2.4) 


The  symbol  A1  will  be  used  to  denote  the  dual  of  A  relative  to  the  bilinear 
form  {  ,  >  ,  that  is 

(Ax,  y)  =  (x,  A'y  >  ,  x  c  a(A)  ,  ye  «(A'). 

The  operator  A'  is  closed  with  domain  jB(A*  )  dense  in  X.  It  is  known 
that  if  A  generates  a  semigroup  S(t)  ,  then  S(t)'  is  also  a  semigroup, 
generated  by  A1  .  See  [4]  for  details  . 

Let  Y  be  a  dense  subspace  of  X  which  is  a  Hilbert  space  in  its 
own  right  with  norm  |  J|y  stronger  than  U  llx  so  ^at  the  Injection 


j  =  Y  — -  X 

1  (y)  -  y ,  yeY, 

is  one-to-one  and  continuous  with  dense  range  Y  C  X  .  We  further  suppose 
that  Y  is  invariant  under  the  action  of  S  (t)'  :  y  e  Y  =>  S(t)'y  e  Y  , 
and  that  this  map  is  continuous  with  respect  to  flS(t)'yUy  ,  ||y  |jy  anc*  1116 

usual  topology  of  [  0,  • )  . 

Let  Y*  be  the  dual  of  Y  with  respect  to  X  as  described,  e.  g.  in 
[1],  [14],  [15].  This  means  that  Y1  is  the  closure  of  X  with 


respect  to  the  norm 


|  (  x,  y  )  | 

JJxJJvi  =  sup  - 

Y  »y«Y 

y/£0  * 


(2.5) 


It  is  known  that  Y1,  so  defined,  is  a  realization  of  the  dual  space  of  Y 
and  it  is  easily  verified  that  the  bilinear  form  (  x,  y  >  may  be  defined,  by 
continuity,  for  x  c  Y1  ,  y  c  Y  as 

<x,  y>  =  lim  ( x.  ,  y ) 
k  —  -  K 

where  {x^}  is  a  sequence  in  X  converging  to  x  in  U  II  y*  *  So 

defined,  (x,y)  generates,  as  x  ranges  over  Y1  ,  all  continuous 

linear  functionals  on  Y  .  We  have 

X  C  X  c  Y1  . 

Definition  2. 1 .  In  the  system  (2. 1)  .  i.e. , 

x  =  Ax+bu,  «  e  L21oc[0,o.  )  , 

b  is  an  admissible  input  element  if  there  exist  Y.  Y1,  as  above,  with 
be  Y1,  such  that  for  every  T  >  0  the  continuous  map 


r 


LT4 


fa 


I*  :  Y  -  C[0,T] 

defined  by 

(Ly)(t)  =  (b,  S  (t)'y  >  ,  yeY,  t  £  [0,  T]  , 
has  a  continuous  extension  to 


(2.6) 


Remark . 
system 


LT  :X  -  LZ[0,T]  . 

It  is  clear  that  this  amounts  to  the  statement  that  in  the  dual  observed 

y  =  Ay 
z  =  (  b,  y  >  , 


b  is  an  admissible  observation  element;  that  is,  for  yeY, 


z(*  )  =  <b,S(.  )'y  >  e  C  [  0,  T  ]  , 

2 

this  relationship  extending  continuously  to  z  (•  )  e  L  [0,T]  for  y  e  X  . 

To  verify  that  Definition  2. 1  enables  consistent  definition,  at  least 
in  a  generalized  sense,  of  solutions  of  (2. 1),  (2.  2)  when  b  is  an 

admissible  input  element  and  to  establish  some  of  the  properties  of  the  resulting 
solution,  we  present 

Theorem  2.2.  If  b  is  an  admissible  input  element,  the  formula 

t  , 

<x(t),y>  =  <xQ,  S(t)'y)  +  /  (b,  S(t-s)  y  >  u(s)ds  ,  yeY,  (2.7) 

0 

defines,  for  each  t  a  0  ,  a  unique  element  x  (t)  e  X  .  Given  T  >  0 
and  u  e  L2[0,  T] 

x(t)  =  S(t)xQ  +  B(t)u  ,  t  e  [ 0, T ]  ,  (2.8) 

where  B(t)  is  the  strongly  continuous  family  of  bounded  operators 
B(t):LZ[0,T]  —X  given  by 

<B(t)u,  y)  =  f  <  b,  S  (t  -  s)'y  >  u(s)ds,  yeY.  (2.9) 

J0 

Proof.  From  (2.  8)  and  the  fact  that  Y  is  dense  in  X  it  is  clear  that 

x(t)  -  S(t)xQ  =  §  (t)  5  B(t)u 

where,  for  yeY, 

<$(t),y)  =  /  <b,S(t-  s)'y  >  u(s)  ds  . 

0 

Let  x  e  X  and  let  {y^}  b®  a  sequence  in  Y  converging  to  x  with 

respect  to  11  *  lx  •  Since  b  is  an  admissible  input  element  the 
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$\ 

l 

i 

i 


i 


corresponding  functions  defined  by 

h^ft-sjs  <b,  S(t-  s)'yk>  (2.10) 

converge  in  L2[0,T]  to  a  function  he  L2[0,T]  .  Defining 

<$  (t),x>  =  /  h  (t-  s  )u(s)ds  , 

0 

we  see  that  for  t  c  [0,T] 

l<e(t>.x>l*|h|  .  u u i  2  , 

L  [0.T]  L*[  0,  T  ] 

*IM  II x  Ux  Mu  II  2, 

I*  [  0,  T  ] 

since  ( cf.  (2.6),  (2.10))  h  =  LjX  .  Hence  ;  (t|  t  X1  =  X ,  This 
also  gives 

lerallx*  ILrBII“ll  2, 

I*  [0,  T] 

showing  that  for  t  c  [  0,T  ]  ,  B  (t)  Is  bounded  with 

||B(t)||  a  ||Lt||  . 

To  establish  that  £  (t)  Is  continuous  In  t  for  each  fixed 
u  €  L  [  0,  T  ]  ( and,  hence,  that  B(t)  is  strongly  continuous  in  t )  ,  let 

A 

OststsT  and  form,  for  y  e  Y 

/N 

<e  <t>-  ?(t),y>  =/  <b,S(t-s)’y)  u(s)ds 

0 

t 

“jf  <b,S(t- s)'y)  u(s)ds  =  (with  x  =  s  -  (t  -  t ) ) 

*  f\b,  S(t- T)y>u(r  +  ^t-t))  d t-/  <b,S(t-s)'y>  u(s)ds 
0  0 
t-t 

+  f  <  b,  S  (t  -  s  )'y  )  u  (s)  ds  = 

0 
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b  L  M 

=  /  <  b,S(t- s)'y)  (u(s  +  (t-t))  -  u(s))ds  +  f  <  b,  S(t- s)'y)u(s)ds 
0  0 

s  liy  II y|lx( ll«(*  +(»-«) -  uii  2  +yu|i  „  ). 

^  X  LZ[0,t]  U  I*  [  o/t  - 1] 

2 

Since  Y  is  dense  in  X  and  since  for  fixed  u  e  L  [0,T]  we  have 

^  Um  Uu||  2  ^  =  0  , 

t-t  —  O'  1T[  0,  t  - 1  j 


linx  U  u  ( -  +(t-t))-uH  ,  =  /v  lim  J|u(»  +  (t-t))  -  u||  ,  =0 

t  t  t  ir[0.t]  t  *  t  ir[ o, t ] 

We  conclude  that  for  fixed  u  e  L  [  0,  T  ]  ,  and  t,  t  as  described , 

unk  iu£>  ‘  i  (t)Uv  =  0 

t  t  t  A  tu 

and  thus  £  (t)  is  continuous  in  X  .  This  completes  the  proof  of  the 
theorem  . 

Let  H  be  a  separable  Hilbert  space  and  let  {pk|  ke  K}  be  a 
sequence  in  H,  K  being  a  countable  ordered  index  set  .  The  p^  are 
strongly  independent  if  no  p^  lies  in  the  closed  span  of  {p^  |  £  /£  k}  . 

If,  in  addition,  there  is  a  positive  number  c  such  that  whenever 

P  =  Z  OfcPfc  »  (2.11) 


the  being  complex  and  Kg  an  arbitrary  finite  subset  of  K  ,  we 


have 


l  i“kr  *°zapiH 

Kn 


(2. 12) 


we  say  that  the  p^  ate  uniformly  I  -independent,  since  (2.12)  implies 

Z  K!2  s  U  p  Ih 


(2. 13) 


.«  w  k  •  r-  •  c 


I 


ig 

I: 

■sS.1 


I 


*;.H 


V>l 

I 


2  v' 

whenever  {a^}  €  *  and  p  =  ^  is  convergent  in  H  . 

K 

If  there  is  a  positive  number  C  such  that 

Ml  =  c2  l  |a.|2 


p  as  in  (2.11),  we  say  that  the  sequence  {p^}  is  uniformly  I  -convergent 
Since  this  property  implies  that  if  {a.}  €  the  series  Va  p,  is 


convergent  in  H  and 

UP#H  *  °2E  l“k|2-  (2-14) 

K 

Recall  that  a  sequence  { p^}  in  H  forms  a  Schauder  basis  for 

H  if  for  every  p  €  H  there  are  unique  coefficients  or ^  such  that  the 

series  E“)cpk  converges  to  p  in  H  ([21]).  A  Schauder  basis 

2 

which  is,  at  the  same  time,  both  uniformly  I  -independent  and  uniformly 
2 

1  -convergent  is  a  Riesz  basis  .  For  evident  reasons  we  shall  also  use, 

2  2 

synonymously,  the  term  uniform  f  -basis.  If  {p^}  is  a  uniform  I  -basis 
for  H  then  every  p  in  H  has  a  unique  convergent  representation 

p  =  £  “kpk 

with  (cf.  (2.13),  (2.14)) 

='2  E I “ki2 *  ll pll h  3  c2Ei«ki2- 

K  K 

For  the  remainder  of  this  section  we  suppose  that 
(i)  the  operator  A  with  dense  domain  fi(A)  c  X  generates  the  strongly 
continuous  semigroup  of  bounded  operators  S(t),  t  »  0  ; 
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(ii)  <J(A)  ,  the  spectrum  of  A  ,  consists  of  discrete,  simple 

eigenvalues  k  e  K  ,  and  the  corresponding  normalized  eigen- 

.  -y 

vectors  k  e  K  ,  form  a  strongly  independent,  uniformly  t 

convergent  Schauder  basis  for  X  . 

Since  the  ^  ,  k  e  K  ,  are  strongly  independent  and  have  closed 
span  equal  to  X  there  exist  unique  biorthogonal  elements  k  c  K  , 

such  that 


<v  v  - 


r  i.  k = i 

\^0,  k  £  l  \ 


k,  f  €  K  . 


As  is  well  known,  the  4^  are  eigenvectors  of  the  dual  operator  A' 
corresponding  to  the  eigenvalues  X.^,  k  €  K  ,  We  further  assume 
(ill)  the  eigenvectors  of  A'  have  the  property 

+k  e  Y  C  X  . 

(this  is  true,  for  example,  if  Y  3  fi((A')r)  for  some  positive  integer  r) 

If  x  e  X  ,  the  fact  that  the  fi  ^  form  a  Schauder  basis  in  X 

implies  the  existence  of  unique  k  e  K  ,  such  that 

K 

the  series  converging  in  X  ,  From  this  it  is  evident  that 

£k  =  x>  »  k  e  K  . 

We  are  not  assured,  in  general,  that  the  §k  are  square  summable  but  the 
uniform  i  -convergence  property  of  the  shows  the  square  summability 

of  the  sequence  {£k>  to  be  a  sufficient  condition  for  convergence  of  (2. 15)  . 


Since  we  assume  the  4*^  lie  in  Y  ,  given  any  element  b  e  Y' 
( and  this  includes  be  X)  we  may  define 

bk  * ( \  • b  > 


(2. 16) 


and  obtain  a  set  of  coefficients  ,  ke  K,  associated  with  b  .  In 
general  it  is  not  possible  to  recover  b  from  the  coefficients  .  (An  exampl 
is  X  s  L2[0,2v]  ,  Y  =  HJ[ 0,  2ir]  ,  «^(x)  =  (2tr )"1eikx  ,  k=  0,1,2,.... 

The  4*jc(  =  here  form  an  orthonormal  basis  for  X  and  belong  to  Y 


but  there  is  a  non- zero  element,  namely  6^  -  5^2tt)  » 


in  Y  for  which 


all  of  the  b^  are  zero .  This  arises,  of  course,  because  the  closed  span 

of  the  4*^  in  Y  is  not  equal  to  Y  . )  As  a  consequence  it  is  not 

generally  meaningful  to  write  b  *  Z  • 

X 

Nevertheless  it  may  be  meaningful  to  consider  the  initial  value  problem 
(2.1),  (2.2),  i.  e. , 

x  =  Ax  +  bu 
x(0)  =  x0eX, 

for  certain  b  e  Y*  ,  namely,  those  that  we  have  already  characterized  as 
admissible  input  elements  .  We  wish  now  to  show  that  the  class  of  such 
admissible  input  elements  can  be  characterized  in  terms  of  the  coefficients  b^ 
and  the  eigenvalues  X  ^  .  If  x  (t)  is  the  solution  of  (2. 1),  (2.  2) 
established  by  Theorem  2.  2  for  an  admissible  input  element  b  then, 
in  particular,  for  t  a  0  , 


<x(t),4^)  *  <xQ,  Sft)1^) 


+  /  (b,  S(t-s)'  4^)  u(s)ds 

V0  rt  Mt-S) 

e  x0  k  +  J  e  u(s)ds 


where 


*0  =  £  *0,^  • 

K  kk' 

We  do  not  know  that  the  numbers  e  x. 


are  square  summable 


but  the  series 


Z8"k\k  K 


must  converge  to  S(t)Xg  by  virtue  of  the  (assumed)  Schauder  basis  property 
of  the  .  It  follows  that  a  sufficient  condition  for  x  (t)  to  belong  to  X 


is  that  the  numbers 


t  Mt-s) 

Ck(t)  =  bk  /  e  K  u  (s) ds 


(2. 17) 


should  be  square  summable  for  each  taO.  Equivalently,  making  a  trivial 
change  of  independent  variable , 

t  Nr8 

^k(t)  =  bkjf  e  f (s) ds  ,  f(s)=  u(t-s). 

The  necessity  of  considering  an  infinite  number  of  values  of  t  can  be  obviated 

2 

by  taking  f  to  be  an  element  of  L[0,T],  T>0  fixed,  and  defining 

f(s)  5  0  in  [t,  T]  for  t  <  T  .  The  map 


rT  V 


Ck=bk  /  8  *(»)<*».  fcl  [0,T], 

*•  o 


so  defined  may  be  designated  as 


^  :  L2[0,T]  -  X, 

LT(f)  =  x  =  S  Ck  ^k  * 
K 


(2. 18) 


(2. 19) 


(2.  20) 


and  it  is  easy  to  see  that  L1^  is  the  dual  of  L,j  :  X  —  L2  [0,T]  as 
defined  by  (2.  6)  .  Thus  the  boundednes  s  of  LT  , 


as  required  in 


Definition  2. 1,  may  be  obtained  as  an  immediate  corollary  if  it  is  shown  that 
L^  ,  defined  by  (2. 18)  -  (2.19)  ,  is  bounded.  For  our  present  purpose 
this  is  the  route  of  choice  . 

Extending  f  further  via  f  (t)  =  0  ,  t  >  T  ,  the  Laplace  transform 
of  f  is  the  entire  function 

OB  T 

j^(z)  =  /  e"ztf(t)dt=  f  e'ztf(t)dt. 

0  0 

In  terms  of  we  clearly  have 

=  bk  ^  *  _X,k^  k  e  K  » 
and  the  following  proposition  is  evident  . 

Proposition  2. 3  .  The  operator  LT  ( equivalently  l't  )  is  bounded 
just  in  case,  for  every  fe  L  [0,T]  the  Laplace  transform  of  f,  has 
the  property 

Elhk^-Vl2  <  -•  (2.  2i) 

K 

We  are  fortunate  that  the  inequality  can  often  be  established  with  the 
use  of  the  concept  of  a  Carleson  measure  and  the  corresponding  Carleson  measure 
theorem  as  it  applies  to  the  space 

=  H^{z  |  Re(z)  >  a)  ,  a  real  .  (2.22) 

2 

The  space  H  {z|Re(z)>  a}  consists  of  those  complex  functions  ^(z), 
analytic  in  Re(z)  >  a  ,  bounded  in  each  half  plane  Re(z)  a  or  +  6,  6  >  0  , 
and  satisfying 

0  *  $  >  a  » 


/  |*M£  +  in)|2dn  =s  M 


(2.  23) 


where  is  a  positive  number  depending  only  on  ^  ( and  not,  in 

particular,  on  |  )  .  It  is  known  (  see,  e.  g.  [10]  )  that  each  such 


function  has  a  limiting  "boundary"  function 


d  (H)  =  lim  +in)  (2.24) 

“  ^  i  a 

defined  almost  everywhere  in  -«»  <  T)  <  ®  and  0))  is  measurable 
with 


f  \taW)\2w  ^  M,  . 

-Q3  ' 


Each  jiS  c  H2 
r  a 

such  that 


is  the  Laplace  transform  of  a  unique  function 


f€  r2loct°-"> 


/  |e'atf(t)|2dt  <  .. 
0 


Let  p.  be  a  (non-negative  valued)  measure  defined  on  the  Borel 
subsets  of  {z  (  z  >  a}  .  Then  p.  is  a  Carleson  measure  if  for  every  real 

t  and  every  h  >  0 


p.  ({  z  j  x -  h  rs  Im  (z)  £  r  +h,  a  <  Re(z)  s  a  +  h})  sAh  (2.  25) 

for  some  positive  A  depending  only  on  p.  ( not  on  h )  . 

For  a  Carleson  measure  we  have 


XiyyjyysJLjL.  If.  P-  is  a  Carleson  measure  on  {z  j  Re(z)  >  a}  with  A 
as  in_  (2.25)  ,  if.  ^  e  ,  and  is  given  by  (2.24),  then 


f  dn  . 

-09 


(2.26) 


i 


A  proof  of  this  theorem  is  offered,  for  the  sake  of  completeness,  in 
Section  4  of  this  paper  .  The  relevance  of  this  theorem  for  our  present  studies 


is  exhibited  in  the  selection  of  a  particular  measure  p  .  For  be  Y'  and  a 
given  discrete  spectrum  {\^}  for  ^  »  *et 

“  *  ’‘Mv,) 

be  defined  by 

=  |bk|2  ,  k  e  K  ,  (2.2?) 

p ({ z |  Re(z)  >  a)  -  { X k  |k  e  K})  =  0  .  (2.  28) 

In  this  case  the  left  hand  side  of  (2.  26)  becomes 

Z  ibk^<-xk}i2  (cf*  (2-21))  * 

K 

The  Plancherel  Theorem,  on  the  other  hand,  gives 

OB  _  00 

/  |^(TD|2dn  =  2ir  /  |e’atf(t)|2  dt 

-w  0 

2S  2e2lalTTT  JT|  f  (t)|2dt 
0 

when  the  support  of  f  is  restricted  to  [0,T]  .  Thus 

Z  2000e2|Q!lT  |/T|f(t)|Zdt 

K  0 

and,  in  view  of  our  earlier  discussion,  we  have 

Corollary  2.  5 .  A  sufficient  condition  in  order  that  b  e  Y*  should  be  an 
admissible  input  element  for  the  system  (2,1)  ,  wherein  cr(A)  =  {X^  |k  €  K} 
and  the  corresponding  eigenvectors  <fi,  ,  k  e  K,  form  a  strongly  independent 


{  z  j  Re  (z)  >  a}  for  some  real  a  . 

We  remark  that  the  assumption  (i)  above  together  with  t 
Theorem  ([4],  [11])  implies  that  the  complex  numbers 
indeed,  confined  to  some  right  half  plane  Re(z)  >  a  .  The 
support  of  f  is  restricted  to  [0,T]  implies  that  the  corresj 
transform  jrf  is  entire  and  satisfies  an  inequality  (2.23)  for 
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Our  first  task  in  this  section  will  be  to  develop  a  method  whereby  input 
elements  b  not  in  the  state  space  X  may  be  identified  as  particular  elements 
of  a  larger  space  Y1  .  The  assumptions  made  will  be  somewhat  more 
restrictive  than  thos  introduced  in  Section  2  .  They  are  by  no  means  necessary 
conditions  . 

Let  us  suppose  that  the  operator  A  ,  generating  a  strongly  continuous 
semigroup  S(t)  on  the  Hilbert  space  X ,  has  (dense)  domain  £>(A)  and 


that  A  possesses  discrete  eigenvalues  ke  K,  with 

lim  |  k|  =  ®. 

p  00  -  • 

Here  p  (k)  denotes  the  number  of  elements  £  e  K  such  that  I  <  k  with 
respect  to  the  assumed  order  relation  on  K  .  The  corresponding  normalized 
eigenvectors  #4^  are  assumed  to  form  a  uniform  basis  for  X  .  We  denote 
the  dual  operator  by  A*  .  It  has  the  same  eigenvalues  and  the  corresponding 

eigenvectors  ^  *  k  e  K  ,  will  be  assumed  normalized  so  that 


<+k-  *i>  * 


1  ,  k  =  £ 

0,  k  /£  £  . 


The  ^  also  form  a  uniform  basis  for  X  ,  as  is  well  known  .  Then  it  is 


easy  to  see  that 


«(A)  =  {y  =  E  xk  **kl  E  I \\ I2  <  “  * 

K  K 

and  that 

fi(A')  =  (y  =  E^k  *kJElxkykl2  <  * 

K  K 

For  the  work  of  this  section  we  take  Y  =  fl(A* )  with  the  graph  norm 


I 


IIyAy  =  Z  i1  +  lxkl2)  iykl2  ' 

K 

where 

CD 

y  =  Z  yk  *k 

k=  1 

in  X  .  Then  Y  C  X  and  the  injection  mapping  is  continuous  .  It  will 
often  be  possible  to  identify  a  Hilbert  space  Z  C  X  with  continuous  injection 
map  such  that  jj  jjg  is  a  familiar  (e.g.  Sobolev)  norm  and  Y  is  a 

closed  subspace  of  Z  on  which  the  norms  j]  jj  z  and  |j  JJ  y  are 
equivalent  . 


We  will  be  concerned  with  two  different  extensions  of  the  operator  A  . 
We  suppose  first  of  all  that  there  is  an  element  x  e  X  not  in  fl(A)  and 
that  L  is  an  operator  on  X  such  that 

<i(L)  =  {£  +ux  |  £  €  £  (A),  u  scalar}  , 

Lx  =  Ax,  x  e  £>(A)  . 


We  will  refer  to  L  as  an  "operational  extension"  of  A  .  Its  significance 
arises  from  the  fact  that  many  of  the  inhomogeneous  boundary  value  problems 
arising  in  applications  can  be  expressed  in  the  form 


dx 

dt 


Lx  , 


(3. 1) 


with  the  restriction 

x  =  £  +ux  €  fi(L)  .  (3.  2 

The  second  extension  of  A  ,  which  is  a  map 

A  :  X  Y*  , 

is  a  standard  one,  often  used,  e.g.  in  [14].  If  y,1)  e  fi(A)  ,  j9(A' )  , 

respectively,  we  have 
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<  Ay,  n  >  =  <  y,  A’  n  >  . 

Since  A*  :  Y  =  fl(A‘ )  ^»  X  is  continuous,  the  form  ( y,  A'tj)  extends 
to  (x,  A'n)  ,  x  e  X  ,  by  continuity  and  density  of  S(A)  in  X  and, 
so  extended,  (  x.  A*  H  )  defines,  for  each  fixed  x  e  X  ,  a  continuous 
linear  functional  on  Y  ,  i.  e. ,  an  element  of  Y1  .  We  define 

A  :  X  —  Y1  =  (  a(A'))' 
by 

<  Ax,  n  >  =  <  x.  A'  T)>  ,  x  e  X,  He  Ys  i(A'), 

Our  first  goal,  with  reference  to  the  system  (3.1),  (3.2),  is  to  replace 

it  by  an  infinite  set  of  scalar  ordinary  differential  equations 


dt 


=  X 


k*k  +  bku'  keK 


(3.3) 


where 


x  (t)  =  Yj  Xjc^)  ^  k  * 
keK 

convergent  in  X  .  In  order  to  do  this  we  recognize  first  of  all  that 

2  *  l  Nc^k 

A 

represents  not  Lx ,  but  rather  Ax  ,  since 

(Ax,  4^.)  =  (x,  A* i|i^>  =  (x,  v|)k>  =  Xkxk  . 

We  rewrite  (3. 1)  in  the  form 

=  Ax  +Lx  -  , 

an  equation  in  y'  .  Then,  since  x  is  to  have  the  form  (3. 2)  with 
£  €  i(A)  ,  and  since 


(3. 4) 


1 


V.v.v.v.  .'.V.  •' 


L£  —  Aj»  =  A£  ,  £  e  fi(A)  , 


(3.4)  becomes 


A  AA 


=  Ax  +  ( Lx  -  Ax )  u  . 


We  define  b  €  Y1 ,  a  continuous  linear  functional  on  Y  =  jB(A'),  by 

...  A  AA  .  .  A  .  A  ■  . 

(b,  T)  )  =  (Lx  -  Ax,  T)  >  =  (Lx,  rj  >  -  (x,  A  tj> 
for  tj  e  jD(A')  =  Y  .  We  then  have 

b  =  ZMk 

K 

where  the  "control  input  coefficients",  bk  ,  are  given  by 

bk  =  (b,i|/k>  =  (loc,  v|/k)  -  (x.A1^) 

=  <Ix,  4^)  -  ^k(  x,  t(/k>  . 

In  most  examples  we  shall  have  iix  =  0  .  Then,  if 

x  =  E  Vk- 

K 

convergent  in  X,  we  obtain,  in  place  of  (3.6), 

bk  =  k  *k  •  k  €  K  • 

Also,  In  this  case,  the  equation  (3.  5)  becomes 

(b,  T  )  =  -(  x,  a'  n> . 

The  equation  (3. 5)  ( or  (3.  8))  will  generally  be  used  to  identify  the 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


functional  form  of  b  while  (3.6)  (or  (3.7))  will  be  used  to  identify  its 
expansion  coefficients  in  terms  of  the  eigenvectors  ^k  of  the  operator  A  . 
While  not  all  admissible  input  elements  can  be  treated  this  way  the  class  is 
large  enough,  we  believe  to  warrant  the  detailed  description  we  have  given  here. 


WWW 


4  •"  •  ^gVjj 


SI 


Example  1 ;  Heat  Equation .  Let  x(s,t)  satisfy 


3x  3 x 

at  n 


0  <  s  <  1  ,  t>0, 


with  boundary  conditions 


x(O.t)  =  0  ,  a  x  (l,t)  +  P  (1,  t)  =  u(t), 
where  <*,0  are  real  numbers,  not  both  equal  to  zero  .  In  this  case 


we  take 


X  =  I*  1 0, 1  ] 

Ax  =  ,  x  €  A  (A)  =  {x  e  IT[0, 1]  I  x(0)  =  0  , 

3s 


Lx  a 


-  &2x 


a x  (1)  +  Px'(l)  =  0}  , 

X  €  0(L)  a  {x  €  H2[0,  1]  |  x(0)  a  0}  , 


A 

X(s)  a 


With 


*cn?r  a  +p  0 

a  +  P  =  0. 


<x,y  >  a  f  x(s)  y  (s)  ds 
0 


we  see  that  if  x,  y  e  0(A) 


(3.9) 


(3. 10) 


(3. 11) 


(Ax,  y>  -  <x,Ay>  af  (x  (s)y(s)  -  x(s)  y"(s))  ds 

0 

1  d 

s  JQ  dl(x'(s)y<s>  '  x(s)y'(s))ds  =  (since  x(0)  =  y(0)  a  0) 
x'(l)y(l)  -  xd)y'(l) 


Ax'd)  +  |x(i)) : 
\x'(l)(y(l) +  §/ 


y(l)  -  X(l)(|y(l)  +y'(l)),  0  /  0 


(1))  -  (x(l)  +|.  x'(l))y'(l),  a  ^  0 


and  we  conclude  A  =  A 


In  the  first  case  of  (3.11)  ,  “  +P  /  0  , 


Lx  s  0  and  we  have,  for  t)  e  fi(A  )  =  (A) 


<b,n>  =  -<x,  A'n)  =  --—p  jf  sn"(s)ds 


,  ,  i  i  -nd)  + 

=  aTp<-sTHs>l0  +  T1<s)l0)  = - at 


-n'(i)  +  n  (i) 
— 


4-  nu),  p  /  o 

p 


-  ffi)'(l)  ,  a  /  o  . 


Thus  we  have 


!5(D’  0 

i  S'  Cl)  .  a  t  o  . 


The  two  agree  if  neither  a  nor  P  are  zero  because  the  linear  functional 
4*  6/i»  +  k  6'(1)  ls  zero  in  (  *>(A))'  =  Y'  in  this  case  • 

P  v1!  a 

2 

The  eigenvalues  of  A  are  where,  for  k  =  1,  2,  3,  . . 

a  sin  ( u  ^)  +  P  « ^  cos  (« ^)  =  0  .  i 


Let 


and  (3.13)  becomes 


J a2 +P2U 

cos (u^- 0^) 


2 

k 


=  sin  0^ 


=  cos  0  k 


0 


so  that 


s 


..VI 


$ 

@1 


“k  =  * — 2 — '  w  +  sin 


Va 2  ♦(»*.*  )  ' 


It  is  easy  to  see  that 


=  6  ( rr)  as  k  — »  •  so 


"k  =  (^T^~  )  w  +  ®  (£  )  .  P  ^  0  , 

“k  =  1  )w+-|skTr»  P  =  0  . 


Defining 


2  r  2 

vk  =  J  sinless)  ds 

it  is  easily  seen  that  in  all  cases  the  are  nonzero  and 


lim  v.  =  — — 
k-«  k  -n/T 


Then  the  eigenfunctions 


2 

form  an  orthonormal  basis  for  L  [0, 1]  .  It  follows  then  that  the  coefficients 
of  the  input  distribution  elements  ( 3. 12)  are  given  by 


"k  = 


sln<V 

w.  cos (u.  ) 


P  A  0. 


a  £  0  . 


(3. 16) 


We  consider  here  the  case  P  /  0  ,  saving  the  analysis  for  3=0 
until  later  in  this  section  .  If  0  £  0  ,  formula  (3. 16)  shows  the  bfc 
to  be  uniformly  bounded  .  The  complex  numbers  -  ^  have  the 

property  (from  (3.14)) 

\  _  i  2k  •  1  |2  2  .  i  <  > 

~  (  2  ^  w  *  ®(1)  • 


(3.17) 


Thus  the  number  of  such  -X^  in  any  set  •  |Im(z)  -  t|  ^  h  , 
asRe(z)s  a  -fh  is  <$  ( and  it  follows  that  the  measure  n  with 
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2 

M-\)  =  |bkr  .  F({Re(z)  a  a>  -  U  {-X^}  )  =  0  is  a  Carleson  measure. 

Hence  if  0^0  the  boundary  input  (3. 10)  is  admissible  . 

In  this  case  the  result  is  easily  obtained  without  the  Carleson  measure 


theorem;  for.  if  the  coefficients 


are  square  summable  and  T  >  0  , 


iZ\%°K\Lz[0'T] 


s*p{\\\}J  l  KI3 

k  k  =  1 


v  rT  2V, 

'  Z  /  e  * 

k=  1  0 


2V 


aj|p{lbkl>Jl  Kl  J  Z  /  8  *  *  s“Ptlb|cl> «/  Z  l°kl 

k—  1  k  —  1  ®  ^  k—  1 


2  TT 

k=i 


(3.18) 

since  sup  { |  bjJ  }  <  •  ,  and  we  conclude  that  the  function  sequence 

{b^e  }  is  f  ^-convergent  in  L^[0,T]  .  Our  next  example  is  chosen 
in  such  a  way  that  a  simple  argument  of  this  type  does  not  apply  and  the  Carleson 
theorem  is  actually  needed  . 


Example  2.  Another  Heat  Conduction  System  .  As  a  further  example 


we  ask  the  reader  to  consider  the  system  shown  in  Figure  3. 1 


Figure  3.1 


The  shaded  horizontal  bar,  B  ,  represents  a  layer  of  material,  whose  depth 
will  be  assumed  negligible,  and  whose  heat  conductivity,  k  ,  is  small  in 
comparison  to  its  specific  heat  R  while  the  region  n  consisting  of  the 
half  strip 

O  :  0  25  X  £  1,  Z  3.  0  , 

is  assumed  filled  with  a  material  whose  specific  heat,  r  ,  is  small  by 
comparison  with  its  conductivity,  K  .  The  heat  flow  equations  are  thus 
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3T  .  a^T  y  a_ 

9t  ax2  32  ' 


(3. 


r  ll  _  K  f  afr  .  _§f_T_  * 
a‘  '  K  ‘  ax2  +  32 2  ’ 

together  with  boundary  conditions 


4l<0't>  =  0,  id.t)  =  o 


U  (0,  z,  t )  =:  0  ,  f£(l,z,t)  =  g(z)u(t), 

(x»  z,t)  =  lim  t  ( x,  z,  t )  —  0 

11m  az  - - * 

t(x,  0,t)  «  T  (x,t)  ,  osxsl  , 


(3. 

(3. 

(3. 

(3. 

(3. 


The  inhomogeneous  boundary  condition  along  x  =  1,  z  :S  0  ,  represents 
the  input  heat  flux  .  In  (3.19),  (3.21),  T(x,  t)  is  the  temperature  in 
the  bar  ,  x(x,  z,t)  the  temperature  in  8  . 

If  we  assume  k,  r  veiy  si.iall  by  comparison  with  R,  K,  we  may, 
as  an  idealization,  replace  f  3-.  19J  and  (J.  20)  by 


8T 

at 


s  »K  — 

*  az 


a2r  ,.32t 

8x2  8z 


r  =  0  ’ 


retaining  the  boundary  conditions  (3. 21)  -  (3.  24) 
state  space 

X  *  {T  *  T  (x)  j  T  €  L2[0,1]  )  . 

We  define  an  operator  A  on  X  with  domain 

«(A)  =  1^(0}  1] 


(3. 

(3. 

We  take  as  our  basic 


WLV.V  JWY 


as  follows  .  Given  T  e  (A)  ,  we  let  t  =  t(x,  z) 

satisfy  (3.  26) 

in  ft  together  with 

t(x,0)  =  T(x),  O^xsl, 

(3. 27) 

and 

0,z)=0,  (1.2)  =  0 

(3.  28) 

lim  •!§(•,  z)=sO  ln  L2[0, 1]  . 

Z  — » 

(3.  29) 

lim  t(-  ,  z)  =  0  in  HX[ 0, 1]  . 

Z  -*  *-  co 

(3.  30) 

From  [14],  for  T  e  H*[0, 1]  we  have  re  H3/2(ft) 

.  The  trace 

theorem  ([1]  ,  [14]  )  then  gives 

•U(.,0)«r  L2[0.1], 

and  we  define 

AT  =s  -  —  —  f  •  0) 

AA  R  3z  1  *  u'  * 

(3.  31) 

So  doing,  (3.25)  becomes 

T*  AT 

(3.  32) 

and  (3. 31)  is  subsumed  in  the  definition  of  A  . 

Lemma  3. 1 .  The  operator  -A  is  the  positive  square  root  of  the  Sturm-Liouville 

S(T)  =  {T  e  H  [  0,  1]  (0)=-^  (1)  =  0}  . 


Proof.  We  compute  (-A)  T  for  T  e  fl(T)  .  For  such  T  the  solution 

of  (3.27)  -  (3.  30)  e  H5/2(ft)  .  If  we  let 


then 


2<x-2>  ■  I  4i  (x-2> 


T  (  %  0)  =  -AT 


2  2 

(-A)2T  =  -At  (•  ,  0)  =  ^  (•  ,  0) 

K  dz * 

...  K2  a2T  ,  K2  d2T 

IT  8x*  dx^ 


since  t  €  H^/2(  ft)  together  with  (3. 26)  implies  that 


2  2 

H(-,0)+-4(',0)S0  in  L2  [  0, 1  ] 
8x  dz 


and  T  =  t  (  •  ,  0 )  . 


The  positivity  of  -A  follows  from  the  divergence  theorem  .  If 
T  c  6(A)  and  t  =  t(x,  z)  is  constructed  as  above,  we  have 

fJU  li(x*z))2  +  (|f  (x,z))2]  dx  dz 

=  //J^7t(x,  z  )||2  dx  dz  (V  =  gradient) 

n 

=  J^f[  dlv  (t(x,z)Vt(x,z))  -  t(x,  z )  a2t(x,  y )]  dx  dz 

(A2  =  Laplacian)  =  (from  (3.26)) 

=  J/div(T(x,  z)  Vt(x,z))  dxdz  =  (using  (3.  27)  -  (3.  30)) 
ft 


This  completes  the  proof  . 


Accordingly,  A  is  self  adjoint  with  eigenfunctions 

j*0(x)  =  1  ,  ^  (x)  =  \Tz  cos  (ktrx ),  k  =  1,  2,  3, . . . 

and  eigenvalues 

X-q  =  0  ,  =  -  ^  k  ir ,  k  =  1,  2,  3, . . .  . 

Let  w(x,  z)  be  the  solution  of  the  following  inhomogeneous 


boundary  value  problem: 


82w  82w 

8x2  8z2 


0  in  Q 


8w 

8x 


(0,z) 


0  , 


lim  -§£-<x,z) 

Z  -<*> 


8w 

8x 


( 1,  z  )  =  g  (z) 


lim  w  { x,  z ) 

Z  -*•  -oo 


0  . 


w  (x,  0 )  =  0  ,  osxsl. 


We  will  assume  that 


g  (z)  is  such  that  the  resulting  w  (x,  z ) 


€  H2(f2) 


In  this  case  the  inhomogeneous  equation  can  be  interpreted  as 
T  =  AT  +bu 

where  b  =  b  (x)  is  given  by 

b  (x)  =  -  ^  ( x,  0 )  . 

To  compute  the  coefficients  of  the  expansion 


b  (x)  =  Yj  \  ^  (x)  » 
k  =0 

we  note  that  since  A  is  self  adjoint,  =  ^(x) , 

bk  =  /  Jrfk(x)  b  (x)  dx  . 


and 


The  Carleson  measure  theorem- can  be  used  in  a  slightly  different  way 

than  that  set  forth  in  Corollary  2.  5  to  show  that  if  g  e  L2  ( -® ,  0  ]  then  the 

2 

are  square  summable  and  b  is,  consequently  an  element  of  L  [0, 1]  , 
Writing  £  =  - z,  g(-£)  =  g(£)»  we  see  that 


bk  = 


-ir  \T  k 

R 


/  exp  ( -kir  £ )  g(£)d£  . 
0 


Since  the  measure  h-  assigning  the  value  1  to  each  of  the  points  kir  , 
k»  0,1,2,,,,,  is  clearly  a  Carleson  measure,  and  since 


changes  only  in  sign,  {b^}  e  I  . 

If  g  (z)  is  Just  bounded  and  measurable  on  -o  <  z  s  0  we  can 
almost  trivially  obtain 


and  the 


bk  = 

will  be  square  summable  . 


It  Is  obviously  possible  to  replace  g  (z)  by  distributions  of  various 


types  .  Taking  g  (z)  =  6, 


corresponds  to  a  point  heat  source  at  the 


comer  x  =  1 ,  z  =  0  and  leads  to 


bk  = 


-1)K  -Jz  K 
R 


(3.  36) 


2  1 

In  our  present  example  X=L[0, 1]  ,  Y  =  $(A)=H[0, 1]  and 

I  * 1 

Y  =  H  [0, 1]  .  The  coefficients  (3,  36)  may  be  recognized  as  those 


corresponding  to  5, 


( referring  now  to  distributions  along  the  x-axis)  , 


Any  measure  p.  assigning  to  the  points  ^  kir  values 

2 

|bjJ  which  are  bounded  evidently  yields  a  Carleson  measure  and  we  conclude 
that  all  of  the  above  cases  correspond  to  admissible  input  elements  .  In  this 
case  the  argument  represented  by  the  inequalities  (3. 18)  will  not  work 

m 

because  the  series  £  (—■ ■  )  is  not  summable  in  this  example  , 


1  2X* 


)  is  not  summable  in  this  example  , 


Example  3  .  Hyperbolic  and  Neutral  Systems 


A  wide  variety  of  systems  involving  linear  hyperbolic  partial  differential 
equations  in  two  independent  variables  x,  t,  or  neutral  functional  equations 
lead  to  systems  of  the  form  described  at  the  beginning  of  this  section,  the 
eigenvectors,  ,  of  A  forming  a  uniform  £  -basis  for  the  state 

space  X  and  the  eigenvalues  X^  confined  to  a  vertical  strip 
at  <  Re  (X)  <  P  in  the  complex  plane.  It  also  usually  turns  out  in  these 
cases  that  the  number  of  X^  in  any  rectangle 

at  <  Re  (X )  <  P,  y  <  Im( X. )  <  6 

is  less  than  or  equal  to  M(6  -  y)  ,  where  M  is  a  fixed  positive  number. 
It  is  evident  that  the  measure  (2.  27),  (2.  28)  is  a  Carleson  measure  in  these 
cases  whenever  the  control  input  coefficients  b^  constitute  a  bounded  set  . 

Example  4  .  Linear  Surface  Waves  .  If  the  operator  A  is  defined  as 
in  (3.  31)  but,  instead  of  the  first  order  system  (3.  32)  we  consider  the  second 
order  counterpart 

i  +  A  £  =  0  (3. 1 

we  obtain  the  linearized  equations  for  small  amplitude  waves  on  the  surface 
of  an  Incompressible  fluid  .  The  theory  is  more  fully  developed  in  [16]  ,  [17] 
[19]  .  With  H  ~  £  ,  (3.  37)  is  equivalent  to  the  first  order  system 


ft)-  ft  •%')■  “ft)  • 


To  obtain  a  topology  corresponding  to  the  energy  of  the  system  one  defines 


-T-%1 


M*  H  =ucr2r  +(Tl.A*TU2r  ,  P.39) 

U\n/lle  L2[0,1]  L2q[0,1] 

where 

Lq [ 0, 1]  =  {£  e  L2[0,1]  |/\(x)dx=  0}  .  P.40) 

0 

2 

The  restriction  to  LQ  [  0, 1]  corresponds  to  conservation  of  fluid  volume. 


On  the  domain 


*0(A)=  k  e  ^tO,  1]  I/1  £(x)dx  =  0} 


the  operator  A  is  invertible  .  Its  eigenvalues  are  (cf.  (3.  34)) 

X^=-^kir,  k=  1,2,3,.., 

with  the  same  eigenfunctions  £$k(x)  *  k  =  1,  2,  3, . . .  ,  as  shown  in 
(3. 33)  .  Correspondingly,  the  operator  G  has  eigenvalues 


(3.41) 


i«k,  -io^,  s  "r" 


J/2  kl/2  _  vk!/2^  k  s  1,  2,  3, . . .  . 


(3.42; 


arid  the  eigenvectors,  orthonormalized  with  respect  to  |[  ||  and  the 

6 

corresponding  inner  product  are 


4*  k  =  (  1  *  k=l,  2,  3, .... 

(3.  43 


To  discuss  admissible  input  elements  in  this  case  we  let  ?k  ,  P_k 
be  non-negative  numbers,  k  =  1,  2,  3, ... ,  and  define 

W-  {i«k)  =  {-i«k}  =  P_k»  k  =  1,  2,  3,  ...  , 

CD 

M.({Re(z)  a  o}  -  U  U  {-i«k}))  =  0  . 


36 


Let  P(w),  -eo  <  a)  <  oo  f  be  defined  as  the  piecewise  linear  function 
such  that  in  the  interval  [  i&j  ^ ,  i  +1  ] 

P(U)  =  PkK-tt"a,)  +  fik-ll<M-“k) 


(3.44) 


"k+1  '  uk 


Since 


r"k+l  1/2 o  \ “lc 2  +  Pk+1  “  k+l 

J**' 2fS(  «)  dcj  =  —-"Y  — 


-  i(Pk  +  *W'  k-»  , 


we  conclude  that  M-  is  a  Carle  son  measure  just  in  case  there  is  a  constant 
C  such  that 

/TcdX/2  P(to)  du  =sC1t-(j|  (3.45) 

a 

whenever  0  <  a  <  t  ,  together  with  a  comparable  condition  involving  the 
and  negative  values  of  u  .  But  (3,  45)  is  true  just  in  case 

"lc  2  pk  S  C  '  k  =  1,  2,  3, . . . 
and  the  comparable  condition  for  negative  k  is 

Uy/  ^  C  ,  ks  1,2,3,...  . 

Thus  for  the  Inhomogeneous  system 


the  input  element 


37 


Is  admissible,  from  this  criterion,  if 

k1/2(|Pk|2  +|P_k|2)  <  C  (3.46) 

for  some  fixed  positive  number  C  .  It  will  be  noted  that  this  is  (slightly) 
less  restrictive  than  the  requirement 


Example  5.  Negative  Results.  For  any  system  similar  to  the  one  in 

-  «k|  =  ®(l/k^^  +  e)  the  Carle  son  measure 
condition  will  be  stronger  than  requiring  b  e  X  .  Hence  failure  of  the 
Carleson  measure  condition  cannot  be  used  to  show  that  an  element  b  is 
not  admissible,  for  any  b  e  X  is  admissible. 

To  illustrate  what  can  be  done  in  a  negative  direction,  we  return  to 
Example  1  with  ^  =  0  .  This  situation  has  been  studied,  using  a  different 
approach,  in  [13]  .  We  present  here  an  argument  more  in  the  spirit  of  the 
present  work.  As  shown  in  (3. 15), 


Example  4  but  with  lwk+i 


and  (cf. 


IT 


2 


(3.16)  and,  w.l.  o.  g.  taking  a  =  1) 


(3.47) 
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b^=^/2kir  cos  (kir  )  =  ( -1  )k  n/2*  kit  .  (3.48) 

Since  Pk  =  lbk|2  =  2k2  tt2  while  (k+1)2  w2  -  k2  irZ  =  2k tt2  +  ir2  , 

2 

it  is  not  hard  to  see  that  the  measure  n  =  =  I  \i  is 

not  a  Carle  son  measure  in  this  case  .  As  we  have  remarked,  this  by  itself 
is  not  enough  to  show  that  the  input  element  with  coefficients  (3.48)  is  not 
admissible  .  To  show  this,  we  ask  the  reader  to  consider  the  function 


k-r 


4ir(z)  =  (z  +1) 

analytic  in  the  complex  plane  minus  the  cut  consisting  of  {z  |  z  real,  z  s  -1}  . 

If  r  >  4- ,  4*  is  square  integrable  on  any  vertical  line  {z  |  Re  (z)  = 

m  r 

£  ,  £  a:  0}  with  uniformly  bounded  L2  norm  and  i| >  r(z)  is  bounded  for 
Re  (z)  a  0  . 


It  follows  that  iji  r(z)  is  the  Laplace  transform  of  a  function 


fr  =  fr(t)  with  f r  €  LZ[0,  «  ) 


Then 


-  k2  2t 

f0  8  V*>  dt 

+  (k2™2 )  = 
r  (k  it  +l)r 

=  M|k|1"Zr),  k  -  «  • 

This  expression  is  not  square  summable  if  r  satisfies  the  inequalities 


(3.49) 


so  we  require 


1  -  2r  i  ^  , 


i  <  r  *  1  • 
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,2  2t 

Let  E  be  the  closed  subspace  spanned  by  the  functions  e  In 

L2[0,«)  and  let  E,j  ,  T  >  0  ,  be  the  subspace  of  L2[0,T] 

A 

consisting  of  restrictions  to  [0,1]  of  functions  in  E  .  If  f  is 
the  orthogonal  projection  of  ff  onto  E  we  clearly  have 

/  e“k  w  *  f  (t)  dt  =  /  e“k  w  *  f  (t)dt. 
r  ^0  r 

It  is  shown  in  [ 7]  ,  [18]  that  the  natural  restriction  map  R:  E  —  E^, 

is  onto,  (obviously)  bounded  and  (not  so  obviously)  boundedly  invertible 

with  respect  to  the  induced  L  [0,  <d),  L  [0,T]  topologies  of  E,  ET 

_k2  2  1 
respectively.  Thus,  with  p^(t)=  e  , 


/V1*2*2*  f(t)dt  *  (t,p  )  2 

0  r  r  kLZ[0,«) 


=  (f  R^Rp.  )  2  =  ((R'Vf.,  RP,.)  2r  , 

r  ^  LZ[0,  <>»)  r  k  LZ[0,T] 


,T  k2w2t 

=  f  e~k  n  1  9  (t)  dt 
J0  r 


where 


9 r  *  (R_1)*tr  e  ET  C  L2[0,T] 


2 

It  follows  that  9f  is  an  element  of  L  [  0,  T  ]  such  that  the  numbers 


.  T  .2  2. 

(-l)K*TZkv  /  e“K  w  r  9  (t)  dt  , 
0  r 


k  s  1,2,3,.. 


are  not  square  summable  .  From  earlier  developments,  the  input  element  b 
with  coefficients  (3. 12)  corresponding  to  the  boundary  condition  (3. 10) , 


with  P  «  0  ,  a  =  1  : 


x(l,t)  =  u(t)  , 


Is  not  an  admissible  input  element  . 
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4.  A  Proof  of  Theorem  2.4.  It  is  clear  that  the  Carleson  measure  theorem 

in  H?,  ,  Theorem  2. 4 ,  is  central  to  our  work  in  this  paper  .  This  result, 

in  one  form  or  another  has  been  known  for  somewhat  more  than  a  decade.  A 
2 

proof  for  H  (D)  ,  where  D  is  the  unit  disc  in  the  complex  plane, 
appears  in  Duren  [  5]  .  A  proof  for  functions  in  is  given,  by 

Koosis  in  his  recent  book  [  12]  .  The  reader  is  also  referred  to  the  recent 
book  [8]  by  J.  Garnett .  Because  the  result  is  not  particularly  well  known 
outside  the  circle  of  mathematicians  working  in  Hp  theory  and  because  the 
results  are  rather  scattered  and  not  readily  available  in  precisely  the  form  we 
require,  we  offer  here  a  proof  of  Theorem  2.4  which  is  a  direct  adaptation 
to  the  half  plane  of  the  result  for  the  unit  disc  appearing  in  Duren' s  book  . 

The  proof  given  here  originally  formed  part  of  the  first  author4  s  doctoral 

\ 

dissertation  [9].  As  in  Duren' s  work,  the  proof  makes  use  of  a  relatively 
simple  case  of  the  Marcinkiewicz  interpolation  theorem  ([23]  ,  Chapter  XI) 
and ,  again  following  Duren,  we  do  not  quote  the  general  Marcinkiewicz 
theorem  but,  rather,  give  a  direct  proof  for  the  simple  special  case  required 
here  . 

We  begin  with  a  covering  lemma  of  "Vitali  type". 

Lemma  4. 1 .  Let  {I^|keA}  =  be  a  family  of  intervals  in  R*  . 

Suppose  there  is  a  positive  number  K  such  that  for  any  finite  collection 

{1^  ,  ,  •••  ,  1^  }  of  disjoint  intervals  in  J 

1  2  n  n 

E  |IKll  <  K.  (4.1) 

k=  1  k 
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Then  we  can  choose  a  sequence  {1^  J  k=  1,2,3,...}  of  disjoint  intervals 
from  3  with  the  property:  for  every  X  €  A  there  exists 
kc  {1,2,3,...}  such  that 


JXC'k 


where  J.  is  the  interval  having  the  same  center  as  Ix  but  five  times 
ic  k 

the  length  of  I.  . 

'  ~  k 

proof.  From  (4.1)  it  follows,  in  particular,  that  the  length,  |IjJ,  of 

1^  Is  uniformly  bounded  (take  n  =  1,  X^  =  X  )  .  Define  the  sequence 

{I.  }  inductively  as  follows  .  Let  I.  be  such  that 

*  k  1 


Ux.l  4.“'!  lJxl  • 


X  e  A 


For  k  *  2,  3.  4, 
and  such  that 


let  I. 


be  disjoint  from  Ix  ,  i=  1,  2,  «.»,n-l. 


|I^  |  a  j;  sup{jl^|  JxeA.I^ni^  =  ^  ,  1*1,  2,  •••,k-l}.  (4.2) 


Since  the  I.  are  disjoint  it  follows  from  (4.1)  that 

\ 


Let  1^  c  J 


llm  |K  |=0. 

k  —  •  k 

Then  there  exists  k  such  that 

n  ixA*. 


(4.  3) 


(4.  4) 


Otherwise  (4.  2)  and  (4. 3)  could  not  both  be  true  .  Let  kQ  be  the 
smallest  integer  such  that  (4.  4)  is  true  .  Then 
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-*rz/ 

Ln=o  |t.T 


|^0(it)|dt 

|t-T|S2N+10 


/  - 

1 1-  t|  ss  a 


M0(it)|dt 


ir?  -l. 

Un=0  2N-; 


?(z)  +  2?(z) 


ProposibLonj4^3.  Let  4*  ^  L^( «  )  and,  for  z  =  o  +  i  t  , 

a  >  o,  let  Iz  be  given  by  (4.5)  while  (cf.  (4.6)) 


$(z)=^sup  yjy  ^  |v|<(t)|dt  . 


(4.8) 


Eet  f1  be  a  Carleson  measure  and;  for  s  2r  0  Jet  E„  be  the  Bore  1 


measurable  subset  of  {z|Re(z)>0}  given  by 
Es  =  (z  |  Re  (z)  >  0  ,  ijT(z)  >  s  >  . 
Then,  with  A  as  in  Definition  2.  3  . 

■kv-Shi.  • 

L  (-»,  00) 


(4.9) 


Proof.  let  j  be  the  family  of  all  finite  intervals  in  R1  such  that 

1  r 

j  1 1  J  l  +  (t)  |  dt  >  s  .  (4.  io) 


^  ^1*  ^2*  *  *  ’  *  *n  6  ^  are  disjoint,  then  (4.10)  gives,  for  every  n. 


n  n 

z  ly z  /  i+widt  hi  j 

=  1  k=l  \  L'f-.-) 


(4.U) 


Thus  j  satisfies  the  hypotheses  of  Lemma  4. 1  and  we  can  find  a  disjoint 
sequence  {In  |  n  s  1,  2,  3, ...  }  c  J  such  that,  Jn  having  the  same 


4»r(t)  = 


0o(it)  if  M0(it)|  >  r 
0  otherwise  . 


2 

From  (i*  )  e  L  (-«,  oo)  ,  we  conclude  that  the  support  of  ipr  is  a 
subset  of  (-00,00)  of  finite  (Lebesgue)  measure  .  Then 

4>r  €  L  (£}r)  C  L  (  £r)  and  we  conclude,  since  ipr  vanishes  outside 
£  ,  that  ipr  e  L^(-oo,oo).  Moreover 

/"ll+rlLl,  dr  =  /"/  |^0(it)|dtdr 


L*(-«  ,=) 


0  Z 


=  /  /  dr|^0(it)|dt  =  /  |^0(lt)|zdt 


‘  I*>“*  >11,2.  , 

L  (“®,  00  ) 


Let  a(s)  =  p(E  )  .  Then  we  can  see  that 

s 

/  ($(z))2d|J.  (z)  =  -/  s2da(s)  =  2 /  sa(s)ds  . 

Re  (z)  >  0  0  0 

From  the  definition  (4.  8)  of  (p  it  is  clear  that  for  any  two  such 

functions,  i|i^,  »  we  have 


Hence 


(+1  +  s  ^(z)  + 
?(z)  =  (+r  +  (^0(i*  )  -+r))(z) 

25  $r(z)  +  (£c*r>  -  Ijf )  (z) 
s  £  (2)  +  r 


since  |^(it)-  »|j  (t)|  is  either  equal  to  0  or  is 

Fs  =  {z  j  £r(z)  >  s}  . 

Suppose  z  e  E2r  .  Then  ££(z)  >  2r  and  (4.15)  gi 

(z)  s  £(z)  -  r  >  r 

and  we  conclude  z  e  F  .  Thus 

r 

E2r  C  Fr  • 

Hence,  from  (4.9)  of  Proposition  4.  3, 

MEr)S|x(Fr/2)S  , 

'  i,  (  -CO,  as  J 


so  that 


/  r«(r  )  dr  =  /  rp(E  )  dr 
0  0  r 

=S  5A  /  |J^  j|  ,  dr  (using  (4.13) 

0  L  (-<*>,  * ) 

*5A|^  (i.JI*, 

L  (-00,  CD  ) 

Then  (4. 14)  gives  the  inequality  (4. 12)  . 

The  proof  of  Theorem  2. 4  is  completed  by  combining 
Proposition  4.  2  with  (4. 12)  above  to  give 

/  |ji*(z)|2dp(z)  =*  ^0  f  (?(z))2d^(z) 

Re(z)>0  w  "Re(z)>0 


* 00 OA  r  i  /  « «i2 
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